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ABSTRACT 

We present the dynamical properties of 97 spectroscopically confirmed ultra-compact dwarfs (UCDs; 
r h A) 10 pc) and 911 globular clusters (GCs) associated with central cD galaxy of the Virgo cluster, 
M87. Our UCDs, of which 89% have M k > 2xl0 6 M 0 and 92% are as blue as the classic blue GCs, 
nearly triple previously confirmed sample of Virgo UCDs, providing by far the best opportunity for 
studying global dynamics of UCDs. We found that (1) UCDs have a surface number density profile 
that is shallower than that of blue GCs in the inner ~ 70 kpc and as steep as that of red GCs at 
larger radii; (2) UCDs exhibit a significantly stronger rotation than GCs, and blue GCs seem to have 
a velocity field that is more consistent with that of the surrounding dwarf ellipticals than with that 
of UCDs; (3) UCDs have an orbital anisotropy profile that is tangentially-biased at radii < 40 kpc 
and radially-biased further out, whereas blue GCs become more tangentially-biased at larger radii 
beyond ~ 40 kpc; (4) GCs with M* > 2xl0 6 M 0 have rotational properties indistinguishable from 
the less massive ones, suggesting that it is the size, instead of mass, that differentiates UCDs from 
GCs as kinematically distinct populations. We conclude that most UCDs in M87 are not consistent 
with being merely the most luminous and extended examples of otherwise normal GCs. The radially- 
biased orbital structure of UCDs at large radii is in general agreement with the “tidally threshed 
dwarf galaxy” scenario. 

Subject headings: galaxies: clusters: individual (Virgo, M87/NGC4486) - galaxies: star clusters: 

general - globular clusters: general galaxies: nuclei galaxies: elliptical and 
lenticular, cD - galaxies: kinematics and dynamics 


1. INTRODUCTION 

Ultra-compact dwarfs (UCDs; Phillipps et al. 2001) 
were originally discovered (Hilker et al. 1999b; Drinkwa¬ 
ter et al. 2000a) as compact stellar systems which are 
more than 1 mag brighter than the known brightest glob¬ 
ular clusters (My ~ —11 mag; Harris 1991) but at least 2 
mag fainter than the prototypical compact elliptical M32 
(My = —16.4 mag). The first five UCDs found in the core 
of the Fornax Cluster were unresolved or marginally re¬ 
solved on ground-based arcsec-resolution images, imply¬ 
ing effective radii of r/,, < 100 pc. Subsequent Hubble 


Space Telescope (HST) imaging of those Fornax UCDs 
(Drinkwater et al. 2004) gave rh > 10 pc, which differs 
significantly from rh of ~ 3 pc for conventional GCs (e.g. 
van den Bergh et al. 1991; Jordan et al. 2005). Since their 
discovery in the Fornax Cluster, similarly bright UCDs 
have been found in other clusters (Virgo: Hasegan et al. 
2005; Jones et al. 2006; Abell S0740: Blakeslee & Barber 
DeGraaff 2008; Coma: Madrid et al. 2010; Chiboucas et 
al. 2011; Centaurus: Mieske et al. 2009; Hydra: Misgeld 
et al. 2011; Antlia: Caso et al. 2013), groups (HCG22 
and HCG90: Da Rocha et al. 2011; NGC1132: Madrid 
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& Donzelli 2013), and even relatively isolated galaxies 
(Sombrero: Han et a. 2009; NGC 4546: Norris & Kan- 
nappan 2011). 

Given the intermediate nature of UCDs, since their 
discovery, there has been ongoing debate about their ori¬ 
gin. The few proposed formation mechanisms in the lit¬ 
erature are: (1) they are merely luminous, genuine GCs 
(Murray 2009), or mergers of young massive star clus¬ 
ters formed in starburst regions, such as those formed 
during collisions between gas-rich galaxies (Fellhauer & 
Kroupa 2002; Bruns et al. 2011; Renaud, Bournaud & 
Due 2014); (2) they are the remains of tidally-stripped 
nucleated galaxies (e.g. Bekki et al. 2003; Goerdt et al. 
2008; Pfeffer & Baumgardt 2013); (3) they are the rem¬ 
nants of primordial compact galaxies (Drinkwater et al. 
2004). 

A consensus about the primary origin of UCDs has yet 
to be reached. In fact, even the name given to this cat¬ 
egory of object, ultra-compact dwarf, has been debated. 
Since the UCD designation implies a galactic origin, they 
have also been referred to as “dwarf-globular transition 
objects” (DGTOs, Ha§egan et al. 2005). In this paper, 
we will refer to these objects as UCDs since that is the 
most common usage in the literature, but our usage is 
not meant to pre-suppose their origin. 

The three properties of UCDs that make them distinct 
from GCs include their larger sizes (e.g. Kissler-Patig, 
Jordan & Bastian 2006), a possible size-luminosity rela¬ 
tion (e.g. Cote et al. 2006; Dabringhausen et al. 2008) 
and slightly elevated dynamical mass-to-light ratios (e.g. 
Ha§egan et al. 2005; Mieske et al. 2008) above a dy¬ 
namical mass of ~ 2xl0 6 M 0 . Mieske, Hilker & Misgeld 
(2012) found that the number counts of UCDs, which 
they defined as stellar systems with My < —10.25, in 
several different environments (the Fornax cluster, Hy¬ 
dra cluster, Centaurus cluster and the Local Group) are 
fully consistent with them being the bright tail of the 
normal GC population. On the other hand, Cote et al. 
(2006) and Brodie et al. (2011) found that UCDs fol¬ 
low dE nuclei, instead of GCs, on the color-magnitude 
diagram, suggesting that most UCDs may be a distinct 
population that is more likely to be related to tidally 
stripped galaxy nuclei, rather than to GCs. 

Recently, Seth et al. (2014) found strong evidence for 
the existence of a supermassive black hole (2.1xl0 7 ) in 
the brightest known UCD - M60-UCD1 (My = —14.2 
mag, Strader et al. 2013), indicating that this UCD is 
most probably a tidally stripped nucleus of a low-mass el¬ 
liptical galaxy. Nevertheless, a spatially resolved analysis 
of the kinematics of the most luminous UCD in the For¬ 
nax cluster (UCD3, My = —13.6 mag; Hilker et al. 1999) 
by Frank et al. (2011) found that its internal kinematics 
are fully consistent with it being merely a massive star 
cluster, without strong evidence for the presence of either 
an extended dark matter halo or a central black hole. 
Moreover, there exists direct evidence that UCD-like ob¬ 
jects can form as supermassive star clusters, such as W3 
in the merger remnant NGC 7252 (M* ~ 7 x 10 7 Mq, age 
~ a few 100 Myr: Maraston et al. 2004) and the recently 
discovered young “UCDs” (Penny et al. 2014) associated 
with star-forming regions in NGC 1275 (a member of the 
Perseus cluster). 

All previous investigations of UCDs were based on ei¬ 


ther incomplete or inhomogeneous small samples, which 
hinders us from understanding the global properties of 
UCDs in any one galaxy or environment. Over the past 
5 years, we have been collecting low-resolution (R ~ 
1300) spectroscopic data for UCDs and luminous GCs 
toward the central regions of the Virgo cluster, using two 
multi-fiber spectrographs: the 2dF/AAOmega (Sharp 
et al. 2006) on the 3.9-m Anglo-Australian Telescope 
(AAT) and Hectospec (Fabricant et al. 2005) on the 6.5- 
m MMT. Our spectroscopic surveys of the Virgo UCDs 
and GCs have been highly efficient (in terms of contami¬ 
nation level of non-Virgo targets), thanks to an unprece¬ 
dentedly clean sample of Virgo UCD and GC candidates 
selected based on the recently completed Next Genera¬ 
tion Virgo Survey (NGVS, Ferrarese et al. 2012), which 
offers deep (c/]i m it = 25.9 mag at lOer for point sources) 
and high resolution (FWHM ~ 0.6" in i band) u*giz 
(and r in the cluster core) imaging data of the Virgo clus¬ 
ter from its core to the virial radius (~ 104 deg 2 ) with 
the MegaCam instrument on the Canada-France-Hawaii 
Telescope. 

Details of the spectroscopic surveys will be presented 
in future papers in a series. In this paper, we present 
a dynamical analysis of the UCDs associated with the 
central cD galaxy M87 {D = 16.5 Mpc; Mei et al. 2007; 
Blakeslee et al. 2009), which hosts the majority of con¬ 
firmed UCDs from our spectroscopic surveys, and thus 
provides the best opportunity for studying the dynamics 
and photometric properties of UCDs as a population. For 
comparison purposes, we also collected radial velocities 
of 911 GCs associated with M87, and did the dynamical 
analysis in parallel with the UCDs. Other papers in the 
NGVS series relevant to the topics covered here include 
a systematic study of photometrically-selected UCDs in 
the three Virgo giant ellipticals M87, M49 and M60 (Liu 
et al. 2015), studies of the distributions of cluster-wide 
GC populations in the Virgo cluster (Durrell et al. 2014), 
a detailed study of the spatial, luminosity and color dis¬ 
tributions of GCs selected based on various NGVS bands 
in the central 2° x 2° around M87 (Lancon et al. in 
preparation), dynamical modeling of M87 GCs (Zhu et 
al. 2014), and the physical classification of stellar and 
galactic sources based on the optical and deep K s imag¬ 
ing (Munoz et al. 2014). Liu et al. (2015) is especially 
complementary to this work, in that it presents a thor¬ 
ough description of the photometry and size measure¬ 
ments of the UCD samples, and a detailed analysis of 
the color-magnitude relation, color distribution, specific 
frequencies, and spatial distribution of the UCDs. 

This paper is structured as follows. In Section [21 we in¬ 
troduce the data and samples used in this work. A brief 
description of the methodology used to select a highly 
clean sample of Virgo UCDs and GCs from the spec¬ 
troscopic catalogs is given in Section [3] The definition 
of our working subsamples is described in Section [4j Sec- 
tion[5]provides an overview of the UCD sample, including 
the spatial distribution, completeness and surface num¬ 
ber density profiles. Section [S] presents the phase-space 
distribution and velocity dispersion profiles of the UCDs 
and GCs. Section Q] presents the velocity distributions 
of our samples. A kinematic modeling of the rotational 
properties of UCDs and GCs is given in Section [H while 
Section [9] is devoted to a Jeans analysis for determining 
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the radial anisotropies of UCDs and GCs. A brief discus¬ 
sion and summary of this paper follow in Section 1101 

2. DATA 

This paper is devoted to a detailed dynamical analysis 
of confirmed UCDs (r^ > 10 pc) associated with the cD 
galaxy M87. To this end, we compiled a sample of spec¬ 
troscopically confirmed UCDs, together with GCs which 
will be used for comparison purpose, from three different 
sources, i.e. our recently finished 2dF/AAOmega AAT 
and Hectospec/MMT surveys, and the radial velocity 
catalogs of Virgo GCs and UCDs compiled by Strader 
et al. (2011, hereafter Sll). For duplicate observations 
among the three sources, a weighted average of the indi¬ 
vidual radial velocities will be used in this work. The lo¬ 
cation of the pointings around M87 covered by our AAT 
(blue dotted circles) and MMT (small red dotted circles ) 
surveys is shown in Figure |T] It can be seen that our sur¬ 
veys covered most of the area encompassed by one scale 
radius (big solid circle ) of the NFW dark matter halo 
toward the Virgo A subcluster (McLaughlin 1999). 

The photometric data for the spectroscopic samples 
are from the NGVS. The reader is referred to Ferrarese 
et al. (2012) for a detailed description of the NGVS. Here 
we only mention that the average seeing of the NGVS g- 
and i-band imaging data in the central 4 square degrees 
around M87 is ~ 0.7" and 0.6", which makes it possible to 
measure (through profile modeling) the size of Virgo ob¬ 
jects down to rh of > 5 pc. The reader is referred to Liu et 
al. (2015, in preparation) for details about size measure¬ 
ments on NGVS images. In addition, we have obtained 
deep Af s -band imaging data for this region (NGVS-IR; 
Munoz et al. 2014). Combining the K s band with the 
NGVS optical band photometry allows us to identify 
most of the foreground stars which would otherwise con¬ 
taminate our spectroscopic catalogs of Virgo UCDs and 
GCs. 

All the magnitudes (in the MagaCam u*griz filters) 
that appear in this paper are on the AB system. In addi¬ 
tion, a subscript of 0 denotes the magnitude in question 
has been corrected for the Galactic extinction. In the re¬ 
mainder of this section, we give a brief introduction to the 
three individual radial velocity catalogs of Virgo UCDs 
and GCs. In addition, we will also introduce the radial 
velocity catalog of early-type dwarf galaxies surround¬ 
ing M87. In this work, the velocity field of surrounding 
early-type dwarf galaxies will be compared to that of the 
UCDs and GCs. 

2.1. The AAT sample 

We have carried out (Mar 28 to Apr 1 in 2012) a sys¬ 
tematic spectroscopic survey of compact stellar systems 
(GCs and UCDs) toward the central regions of the Virgo 
cluster (Virgo A subcluster), using the 2dF/AAOmega 
multi-fiber spectrograph on the 3.9-rn Anglo-Australian 
Telescope (AAT). The survey consists of nine 2dF point¬ 
ings, covering a total sky area of ~ 30 deg 2 . The 8 point¬ 
ings around M87 are shown as blue dotted circles in Fig¬ 
ure [TJ The observations cover a wavelength range from 
~ 3700 A to 8800 A, with a resolution of R = 1300. 

The candidates for our survey were selected to fall in 
the region occupied by spectroscopically confirmed Virgo 
GCs and UCDs on the MegaCam u* — g vs. g — i dia- 



Fig. 1.— Location of the fields covered by our AAT ( blue dotted 
circles , 1° in radius each) and MMT ( small red dotted circles , 0?5 
in radius each) surveys. The big black solid circle marks the scale 
radius (2? 143 = 0.617 Mpc) of the NFW model fitted to the sur¬ 
rounding dark matter halo of the Virgo A subcluster (McLaughlin 
1999). The black plus marks the photometric center of M87. 

granfl, and have 18.5 < g < 20.5 magnitude (—12.6 < 
M g < —10.6 at a distance of 16.5 Mpc for the Virgo clus¬ 
ter). The Virgo UCD candidates were mainly selected to 
have 10 < rh.NGVS ^5 30 pc, and the GC candidates have 
1 < Di.ngvs < 10 pc, where rh.NGVS is the half-light ra¬ 
dius (assuming a distance of the Virgo cluster) measured 
based on the NGVS g and i images. 

The AAT survey obtained radial velocities of 55 Virgo 
UCDs and 52 GCs, of which 22 UCDs and 20 GCs have 
no published velocities before. At the limiting magnitude 
of g ~ 20.5, we obtained radial velocities with errors of 
~ 30 km s -1 in a typical exposure time of 1.5 hour. 

2.2. The MMT sample 

In 2009 and 2010, we used the Hectospec multifiber 
spectrograph (Fabricant et al. 2005) on the 6.5-meter 
MMT telescope to carry out an extensive spectroscopic 
survey (3650-9200A, R = 1000) of the central 2° x 2° 
(576 x 576 kpc) around M87 in three observing runs. 
The MMT pointings around M87 are shown as small red 
dotted circles in Figure [I] 

Similar to the AAT survey, the GC and UCD candi¬ 
dates for this survey were selected using MegaCam u*giz 
photometry from NGVS imaging. At the limiting magni¬ 
tude (g < 22.5) of this survey, we obtained radial veloc¬ 
ities with errors of ~ 30 km s -1 in two hour exposures. 
This survey produced radial velocities for 324 GCs and 
51 UCDs, of which 207 GCs and 18 UCDs (excluding the 
ones discovered by our AAT survey) have no published 
velocities before. 

* By the time we prepared input catalogs for our spectroscopic 
surveys (including the MMT survey described below), the NGVS- 
IR K s band data, which is very efficient in separating out the 
foreground stars, was not available. 
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2.3. The Sll sample 

By combining radial velocities from the literature 
(Huchra & Brodie 1987; Mould et al. 1990; Cohen 2000; 
Hanes et al. 2001; Jones et al. 2006; Hasegan et al. 2007; 
Evstigneeva et al. 2007; Firth et al. 2009; Paudel et al. 
2010) and from their new observations, Sll compiled a 
sample of 927 radial velocities toward the central ~40' of 
M87. Of these, 737 were classified as Virgo GCs and 
UCDs. Given the high-quality multi-band imaging data 
from the NGVS, we re-classified (Section [3]) the original 
927 objects from Sll, and found another 5 low-velocity 
objects belonging to the Virgo cluster. 

2.4. Earty-type Dwarf Galaxies Surrounding M87 

Our surveys extend to the Virgo intra-cluster region. 
We will compare the kinematics of UCDs and GCs in 
the outermost part of M87 to the surrounding early-type 
dwarf galaxies, in order to explore any possible connec¬ 
tion. Within a geometric radius of 2° from M87, there 
are 326 galaxies (either with radial velocity unavailable 
or < 3500 km s -1 ) classified as either dE or dSO galaxies 
in the Virgo Cluster Catalog (VCC, Binggeli & Cameron 
1991). Among the 326 galaxies, 67 were further classified 
as nucleated dE galaxies (dE, N)0. 59 of the 67 dE,N 
galaxies and 67 of the non-nucleated dE/dSO galaxies 
have radial velocities available in the literature, as com¬ 
piled by the GOLDMine project (Gavazzi et al. 2003). In 
addition, our AAT and MMT surveys obtained the first 
radial velocity measurements for another 2 dE galaxies, 
i.e. VCC1317 (V = 327±39 km s” 1 ) and VCC1244 (V = 
824±33 km s -1 ). In this work, we will be comparing the 
velocity field (number density profiles) of the 128 (326) 
galaxies with that of the M87 UCDs and GCs. 

3. IDENTIFICATION OF VIRGO OBJECTS 

3.1. Culling Out the Virgo Objects 

There exists contamination from both background 
galaxies and foreground stars in our spectroscopic cata¬ 
logs. There is a well-defined gap between the Virgo galax¬ 
ies and background galaxies at radial velocities ~ 3000 
km s , so a simple cut in radial velocity at 3000 km s^ 1 
can remove all background galaxies. At the low-velocity 
end (Uos < 400 km s -1 ), to cull foreground stars from 
the spectroscopic catalogs, we made use of the u* — i 
vs. i — K s color-color diagram, which has been shown 
to clearly separate nearly all foreground stars (Munoz et 
al. 2014) from Virgo stellar systems (with the exception 
of some metal-poor G-type stars). For sources that fall 
inside the overlap area of the bona-fide Virgo members 
and foreground stars on the u* — i vs. i — K s diagram (see 
Munoz et al. 2014), we further require that the sources 
should have half-light radii rh measured (on the NGVS 
g- and /-band images) to be > 0.06", corresponding to 
a linear scale of ~ 5 pc at the distance of Virgo clus¬ 
ter. The high-quality NGVS g and (especially) i imaging 
data can resolve Virgo sources down to rh ~ 5 - 10 pc. 

The half-light radius rh,NGVS of each source was de¬ 
rived as a weighted average of two independent measure¬ 
ments in the g and i bands by fitting PSF-convolved King 

t One should keep in mind that the real fraction of dE, N galaxies 
is most probably much higher than 21% (67/326), as demonstrated 
by Cote et al. (2006) based on high-quality HST imaging data of 
100 early-type galaxies in the Virgo cluster. 


(1966) models to NGVS images with the KINGPHOT 
software package (Jordan et al. 2005). Briefly, when us¬ 
ing KINGPHOT, we adopted a fixed concentration pa¬ 
rameter c of 1.5 and a fixed fitting radius r^t of 1.3" (~ 
105 pc at the Virgo distance). According to Jordan et al. 
(2005), the KINGPHOT size measurement suffers from 
large biases when rh > rg t / 2, which is however not ex¬ 
pected to be a problem for our analysis because all but 
one (VUCD7) of previously confirmed Virgo UCDs have 
r h < 50 pc. 

3.2. Separating the Virgo UCDs from GCs 

A UCD is defined to have 10 < rh < 100 pc in this 
work (see Section [4]). At 10 < rh < 20 pc, the Virgo 
UCDs are only marginally resolved in the NGVS images. 
Assuming a Gaussian-shaped PSF, the typical FWHM of 
the NGVS /-band seeing disc (0.6") is equivalent to an rh 
of 0.172", which corresponds to ~ 14 pc at the distance 
of the Virgo cluster. Therefore, size measurements based 
on the NGVS images are especially sensitive to the S/N 
and possible inaccuracy of the PSF, and are unavoidably 
subject to relatively large uncertainties compared to the 
measurements based on HST images. Our test with the 
NGVS images suggests that, sources with g > 21.5 - 22 
mag are subject to relatively large bias and uncertain¬ 
ties (> 20%) in their size measurement, and thus are not 
suitable for our analysis. To pick out a clean sample of 
UCDs based on NGVS images, we require UCDs to have 
Hi, ngvs > 11 pc, Ar h ,NGVs/ni,NGVS < 0.1, and g < 21.5 
{M g < —9.6). In addition, sources with rh.HST > 9.5 pc 
based on measurements with existing HST images (e.g. 
Sll; Jordan et al. 2005) are also included as UCDs, re¬ 
gardless of their brightness. All the other confirmed Virgo 
compact clusters are regarded to be GCs. 

By comparing our size measurements with that deter¬ 
mined with existing HST imaging data (see Table 0, 
our size criteria of UCD selection based on the NGVS 
measurements result in zero contamination from Virgo 
objects with rh.HST < 10 pc. Among the old sample of 
34 Virgo UCDs with c^hst > 9.5 pc, 3 did not have 
NGVS size measurements due to their proximity (within 
10") to saturated foreground stars, 28 have n^NGVS > 
11 pc, 1 (T15886: g = 22.97 mag) has 10 <rh ngvs< 11 
pc, and the other 2 (S6004: g = 21.32; S8006: g = 20.53) 
have rh measured to be less than 10 pc either in NGVS 
g or i band. Therefore, by selecting UCD-sized objects 
( r h > 10 pc) based on the NGVS images, we may miss 
~ 6% of genuine UCDs at g < 21.5 mag and < 3% at g 
< 20.5 mag. 

3.3. The Samples of UCDs and GCs 

Given the above selection procedure, we end up with a 
total number of 97 UCDs and 911 GCs, which fall within 
1.5° of M87 and are not associated with any galaxies 
other than M87 based on spatial location and radial ve¬ 
locities. Some of these UCDs and GCs, especially those at 
the outermost radii, probably belong to the intra-cluster 
population. The full sample of UCDs is listed in Table 
0 The sample of GCs has been recently used by Zhu 
et al. (2014) to determine the dynamical mass profile of 
M87. The full catalog of GCs will be presented elsewhere 
(Peng et al. 2015, in preparation). The 34 UCDs con¬ 
firmed previously in the literature (old sample, Brodie et 
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al. 2011) and the 63 newly confirmed UCDs (new sam¬ 
ple) are listed separately in Table [5] For Virgo objects 
that were already spectroscopically confirmed (as com¬ 
piled by SI 1), we follow the old naming; for the newly 
confirmed Virgo members, we adopted a naming scheme 
which starts with “M87UCD-”. Note that, the column 
r h,NGVS gives the weighted average half-light radius mea¬ 
surements in the NGVS g and i bands. The t^hst mea¬ 
surement (Jordan et al. 2005; Jordan et al. 2009; SI 1), if 
available, is also listed. We point out that the uncertain¬ 
ties of rh.NGVS reported in Table [5] only include the for¬ 
mal errors returned from the KINGPHOT fitting, and do 
not take into account any potential systemic uncertainty, 
such as the degree of accuracy of the PSF and suitability 
of King models for representing the UCD light profiles. 

Our sample of UCDs is nearly 3 times larger than pre¬ 
viously known (34) Virgo UCDs, and our GC sample is 
~ 20% larger than that of Sll. Radial velocities for 39 of 
the 97 UCDs were obtained from our AAT and MMT sur¬ 
veys for the first time. In terms of spatial coverage, one 
of the most important improvement of our GC sample is 
at the projected galactocentric distances larger than 30' 
from M87. Specifically, our sample includes 63 GCs in 
the projected radius range from 30' to 60', and this is 7 
times larger than that of Sll. While the full catalog of 
M87 GCs will be presented elsewhere, we emphasize that 
the main results related to GCs in this paper would not 
change qualitatively if only the Sll sample of GCs was 
used in our analysis because of the already large spec¬ 
troscopic sample of M87 GCs in the literature. Since the 
surface number density of UCDs is relatively low, our 
kinematic analysis will be carried out in a coarser spatial 
(or radial) resolution than previous studies (Cote et al. 
2001; Sll). Therefore, whenever relevant, we refer the 
readers to Cote et al. (2001) and Sll for a more detailed 
kinematical analysis of M87 GCs within the central 30'. 

4. DEFINITION OF WORKING SAMPLES 

The primary goal of this work is to explore the dif¬ 
ferences or similarities between UCDs and GCs. To this 
end, we define the following subsamples in this paper. 

• UCDs. UCDs are distinguished from GCs as hav¬ 
ing half-light radius rh > 10 pc (although in prac¬ 
tice we require rh,NGVS > 11 pc, as described 
above) in this work. The two most commonly 
adopted definitions of UCDs are mass (2xl0 6 < 
Mdyn < 10 8 M 0 ; e.g. Hasegan et al. 2005; Mieske 
et al. 2008) and/or r h (10<rh<100 pc; e.g. Norris 
et al. 2011; Brodie et al. 2011). The mass defini¬ 
tion is justified by the findings that 1) compact 
stellar systems with M > 2x 10 6 tend to have M/L 
significantly higher than the lower mass systems; 
2) there seems to be a size-luminosity relation set¬ 
ting in at M > 2xl0 6 (e.g. Rejkuba et al. 2007; 
Evstigneeva et al. 2008; Dabringhausen et al. 2008; 
Norris & Kannappan 2011), in contrast to the more 
or less constant rh (~ 3 pc, e.g. van den Bergh et 
al. 1991; Jordan et al. 2005) of “normal” GCs. The 
size definition of UCDs differentiates them from 
normal GCs as dynamically un-relaxed stellar sys¬ 
tems (e.g. Mieske et al. 2008). The two definitions 
may converge at the highest mass end. While we 
adopted the size definition in this work, we will try 


to explore the significance of mass in differentiating 
UCDs as stellar systems distinct from normal GCs. 

• Blue GCs and Red GCs. A double-Gaussian 
fitting to the NGVS (g — i) o bimodal color dis¬ 
tribution of photometrically-selected GCs in M87 
suggests that the blue and red components cross at 
( g — i) o = 0.89 mag. Therefore, we classified GCs 
as blue (V=683) and red (N= 228) at a dividing 
(g — i) o= 0.89 mag. 

• Bright GCs and faint GCs. Bright GCs and 
faint GCs are separated at NGVS io = 20.5 mag¬ 
nitude, which corresponds to a stellar mass of ~ 
1.6 - 2x10 6 Mq at [Fe/H] ranging from —1.3 (the 
typical value for blue GCs in M87; Peng et al. 
2006) to —0.3 (the typical value for red GCs in 
M87) for a 10 Gyr old stellar population with a 
Chabrier or Kroupa stellar initial mass function 
(IMF). The dividing magnitude (or mass) was cho¬ 
sen to roughly correspond to the proposed mass 
boundary between UCDs and GCs for the mass 
definition of UCDs. By separating the bright GCs 
from the faint GCs, we will explore the importance 
of mass or luminosity in differentiating UCDs from 
GCs. 

The median (g — i) o of our samples of UCDs and blue 
GCs are 0.75 and 0.74 respectively, and about 92% of our 
UCDs fall into the color range of the blue GCs. So we 
will place additional emphasis on a comparison between 
dynamical properties of UCDs and blue GCs throughout 
this paper. 

5. OVERVIEW OF THE UCD AND GC SAMPLES 
5.1. 2D Spatial Distribution 

Figure [2] presents the spatial distributions of all spec¬ 
troscopically confirmed UCDs, GCs and dE galaxies 
around M87. Data from our AAT survey, MMT sur¬ 
vey, and Sll are represented by different colors. Different 
types of objects are plotted as different symbols. When 
plotting Figure [2] for duplicated observations among the 
three data sources, we group them into, in order of pri¬ 
ority, the Sll catalog, the AAT survey catalog, and the 
MMT survey catalog. Figure [3] shows the spatial distribu¬ 
tion of UCDs, color-coded according to their line-of-sight 
velocities. 

Like the GC system (e.g. McLaughlin et al. 1994; Forte 
et al. 2012; Durrell et al. 2014), the spatial distribution of 
the UCDs broadly follows the stellar diffuse light. Within 
2° of M87, the outermost confirmed UCD has a projected 
galactocentric distance R v of 55' from the center of M87, 
the outermost red and blue GCs have R p = 49' and 
85' respectively, with only one GC lying within 65' < R p 
< 85'. 

An interesting “overdensity” of ~ 11 UCDs can be seen 
toward the northwest of M87 between ~ 15' and 30'. 
After checking the radial velocity distribution of UCDs 
belonging to this spatial “overdensity”, we found that 
the “members” of this overdensity have radial veloci¬ 
ties ranging from ~ 900 to 1750 km s -1 , suggesting this 
“overdensity” is due to chance alignment, rather than a 
physical substructure. 
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• AAT GCs 

Sll GCs 

• MMT GCs 
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o AAT UCDs 
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Fig. 2.— Spatial distribution of spectroscopically confirmed GCs, UCDs, and dE galaxies within the central 1.5° (left panel) and the 
central 20' (right panel) of M87. Dotted ellipses represent the stellar isophotes of M87 at 3 R e , 5 R e , 10 R e and 20-R e3 where R e is the effective 
radius measured based on the NGVS g-band image. Different colors indicate data sets from different sources. The large black plus marks 
the photometric center of M87. 
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Fig. 3.— Spatial distribution of spectroscopically confirmed UCDs around M87. As in Figure |2j dotted ellipses represent the stellar 
isophotes of M87 at 3 R e , 5 R e , 10 R e and 20R e . The plotted objects are color-coded based on their line-of-sight velocities. 
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In what follows, we will be mostly working with the 
geometric average radius R av when exploring various ra¬ 
dial trends. R av is defined to be equal to a \/1 — e, where 
a is the length along the semi-major axis (PA ~ 155°) 
and e is the ellipticity. Given that the spatial distribution 
of UCDs and GCs, in terms of flattening and orientation, 
roughly matches the stellar diffuse light of M87 (Durrell 
et al. 2014), we adopted the radial profiles of e and PA 
of the stellar isophotes of M87 determined with the high- 
quality NGVS g-band imaging data. Our measurements 
of e, ranging from ~ 0 in the central 0.5' to 0.33 around 
10' 15' along the semi-major axis, are in good agreement 
with previous studies (e.g. Ferrarese et al. 2006; Kor- 
rnendy et al. 2009). Measurements of e beyond ~ 15' are 
subject to relatively large uncertainties, so we fixed e as 
0.33 at a > 15'. 
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Fig. 4.— Completeness levels of our sample of UCDs as a func¬ 
tion of magnitude in g (black solid lines) and i (red dotted lines) 
bands. 
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5.2. Completeness of the UCD Sample 

To quantify the completeness of our spectroscopic sam¬ 
ple of UCDs, we selected all of the UCD candidates 
within the central 60' (in geometric average radius) of 
M87 from the NGVS photometric catalog (Munoz et al. 
2014), based on the u* — i vs. i — K s diagram (Section[3|), 
size measurement (ll<?'h,NGVS <50 pc), and isophotal 
shapes (SExtractor ellipticities < 0.25). In addition, 
any obviously extended galaxies were further removed 
from our final photometric sample. 

The completeness in a given magnitude range was de¬ 
termined as the j^ spe ° x 100, where N can o is the number 
of photometric candidates and N spec is the number of 
candidates with radial velocity measurements. The com¬ 
pleteness in different go and i o magnitude bins is shown 
in Figure |U Overall, our sample is expected to be ~ 60% 
complete at go < 21.5 mag, and ~ 55% complete at io < 
20.5 mag. In particular, the sample is ~ 98% complete 
at go < 20.5, which corresponds to M g < —10.6. 


TABLE 1 

Sersic Profile Fitting 


ID 

(1) 

90 

(mag) 

(2) 

N e 

(arcmin -2 ) 

(3) 

Re 

(arcmin) 

(4) 

n 

(5) 

UCDs 

18.5 - 20.5 

0.03 

12.57 

1.43 

Blue GCs 

18.5 - 24.0 

0.87 

20.10 

3.03 

Red GCs 

18.5 - 24.0 

3.61 

5.89 

4.26 


Note. — (1) Population name. (2) g-band magnitude 

range of the sample. (3) Surface number density at the effec¬ 
tive radius R e . (4) Effective radius. (5) Sersic index. 


Figure [5] presents the io magnitude distribution as a 
function of R av from M87. The UCDs, blue GCs and 
red GCs are plotted as different symbols in Figure [5] We 
can see that the available observations of GCs within the 
central ~ 5' reach down to 22 mag, which is about 1 mag 
fainter than that at larger radii. We point out that ~ 89% 
of our UCDs have io < 20.5 mag. 

5.4. Surface Number Density Profiles 

As is shown in Figure |H our sample of UCDs is ex¬ 
pected to be nearly 100% complete down to go < 20.5 
mag. In Figure[6l we show the radial number density pro¬ 
file of the 59 UCDs with g 0 < 20.5 mag, together with 
profiles of the blue GCs and red GCs (18.5 < go < 24.0 
mag) determined by Durrell et al. (2014). Number den¬ 
sity profiles for the surrounding dE galaxies were also 
shown for comparison. Our specific choice of radial bin¬ 
ning for constructing the profile of UCDs ensures that at 
least 10 data points fall into each radial bin. The verti¬ 
cal error bars of the radial profiles represent the poisson 
noise. Note that radial profiles of the GCs have been ver¬ 
tically shifted down arbitrarily (2.1 for the blue GCs and 
1.7 for the red GCs) for comparison purposes. We point 
out that the GC surface density profiles derived by Dur¬ 
rell et al. (2014) were based on an adaptive-smoothed 
GC density maps, with the smoothing kernel FWHM > 
3' - 5'. Therefore, the intrinsic surface profiles of GCs in 
the inner radii may be a little steeper than those shown 
in Figure [6] 

We adopt the Sersic function (Sersic 1968; Ciotti 1991; 
Caon et al. 1993; Graham & Driver 2005) to quan¬ 
tify the radial profiles. Due to the small sample size of 
sparsely-distributed UCDs, instead of simply fitting the 
binned profiles, we used the maximum likelihood method 
(e.g. Kleyna et al. 1998; Westfall et al. 2006; Martin et 
al. 2008) to estimate the UCD number density profile. 
Specifically, the likelihood function to be maximized is 
defined as 


£(N e ,R e ,n) oc Y[ii(Ri\N e ,R e ,n) (1) 

i 

where £i(ri\ _/V e , R e , n) is the probability of finding the 
datum i at radius Ri given the three Sersic parameters, 
i.e. the effective radius R e , the number density N e at R e , 
and the Sersic index n. In particular, 


5.3. Magnitude vs. Galactocentric Distance 


fi {R'i | N (., R v . 77.) oc N e exp \ h Tl 




( 2 ) 
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Fig. 5. — Radial distribution of iq magnitude for UCDs and GCs. 
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Fig. 6.— Radial surface number density profiles M87 UCDs 
(black) with go < 20.5 mag. The radial profiles of blue GCs (blue) 
and red GCs (red) determined by Durrell et al. (2014) were also 
plotted for comparison. Also over plotted are the number density 
profiles of surrounding dE galaxies. The profiles of blue GCs and 
red GCs have been vertically shifted arbitrarily for comparison 
purpose. Overplotted on the data are the best-fit Sersic profiles for 
UCDs, blue GCs and red GCs. Note that our sample of UCDs with 
go < 20.5 mag is expected to be nearly 100% complete. 


where b n is a constant that is defined as a function of n 
such that R e is the effective radius, and we adopted the 
formula determined by Ciotti & Bertin (1999) to relate 
b n to n. Furthermore, N e can be expressed as a function 
of R e , n and the total number of UCDs N tot (= 59) 
by integrating the Sersic profile over the projected area 
nR 2 to the limit radius R\i m = 60' (e.g. Graham & Driver 
2005). Specifically, 




N t< 


2TrR%n ex ^ n ^ j(2n, b n (Ru m /Re) 1/n ) 


(3) 


where 7(2n, b n (R\\ m /R e ) l t n ) is the incomplete gamma 
function. By substituting N e from Equation [3] in Equa¬ 
tion [1] £ was maximized to find the most likely parame¬ 
ters R e and n (and thus N e ). For the blue and red GCs, 
we directly fit the Sersic function to the binned radial 
profiles shown in Figure [H which is adequate given the 
large sample size of the photometric samples of GCs. 

The best-fit Sersic profiles of the three populations are 
overplotted in Figure[6l and the most likely estimation of 
the Sersic parameters is listed in Table [I] The difference 
between radial profiles of UCDs and GCs is significant. 
The UCDs have the shallowest radial profiles in the in¬ 
ner ~ 15' among the three populations, and in the outer 
radii the profile of UCDs is as steep as that of the red 
GCs. Previous studies have shown that the surface num¬ 
ber density profile of the red GCs closely follow that of 
the diffuse stellar light (e.g. Geisler et al. 1996; Harris 
2009; Durrell et al. 2014). In addition, the surrounding 
dE galaxies have much flatter and extended number den¬ 
sity profiles than UCDs and GCs. 

6. PHASE-SPACE DISTRIBUTION AND Vrms PROFILES 

Figure [7] gives the radial variation of the line-of-siglit 
velocities of UCDs, GCs and dE galaxies. The blue and 
red GCs are plotted separately. A ^-magnitude color- 
coded plot of Vos vs. i? av for UCDs is shown in Figure 

m 

While we will present a detailed kinematic modeling of 
the rotation and intrinsic velocity dispersion for our sam¬ 
ples in Section [8l it is helpful to first explore variations 
in the root-mean-square line-of-sight velocity u rms as a 
function of the galactocentric distance from M87 (Figure 
®. To construct the u rms profiles, we adopted the “slid¬ 
ing bin” method. Specifically, bins of fixed radial width 
were slid from the center outward, with an offset of 1' be¬ 
tween adjacent bins. u rms was calculated for data points 
falling into each individual sliding bin. Considering the 
lower number density of UCDs and GCs in the larger 
radii we used different bin widths in different ranges of 
radii. The bin widths are: A R av = 4' for R av < 8', A R av 
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= 8' for 8'< i? av < 20', A R av = 12' for 20'< i? av < 40', 
and Ai? av = 16' for R av > 40'. 

We require at least 19 data points to be available in 
each bin when constructing the profiles shown in Figure 
[9] A biweight v rms was calculated for each sliding bin 
following the methodology of Beers, Flynn, & Gebliardt 
(1990; cf. Equation 9). The biweight ty ms estimator is 
relatively outlier-insensitive, and has proven to be supe¬ 
rior to the classical formula when dealing small samples. 
The 68% confidence interval for v lms was estimated by 
randomly resampling the real data sets without replace¬ 
ment. t; rms in some fixed radial bins is also overplotted in 
Figure [21 The dispersion profile for the surrounding dE 
galaxies is also shown for comparison. 

As shown in Figure [9] the UCDs follow a velocity dis¬ 
persion profile more similar to that of the blue GCs than 
the red GCs. The i; rms of blue GCs beyond 30' increases 
steeply to reach ~ 500 km s -1 . The red GCs have an 
overall lower velocity dispersion than the UCDs and blue 
GCs. The blue GCs and UCDs show a slight increase in 
velocity dispersion between ~ 4' and 12', whereas the 
red GCs do not clearly show such a “hot” feature. This 
may suggest the existence of “hot” substructures which 
have not yet reached an dynamical equilibrium state, as 
suspected by Zhu et al. (2014). 

We note that the rising u rms of the blue GCs beyond ~ 
30' should not be regarded as signifying a transition to 
the general cluster potential, because the rising disper¬ 
sion is mostly driven by a larger scatter of Vios toward 
the low-velocity side of the systemic velocity of M87 (see 
Figures [7] and [TT]), and the majority of the blue GCs out 
to R ~ 40' (~ 190 kpc) are actually tightly clustered 
around the systemic velocity of M87. A low-velocity ex¬ 
cess (toward the North West of M87; Figure H9jl is also 
present in the velocity distribution of dwarf galaxies lo¬ 
cally surrounding M87 (Binggeli, Popescu & Tammann 
1993). So an excess of low-velocity GCs beyond ~ 30' in¬ 
dicates a significant contamination from the intra-cluster 
population of GCs projected along the line of sight. 

Moreover, beyond the central ~ 30', the GCs that are 
clustered around the systemic velocity of M87 tend to 
have redder colors than those with lower velocities. If 
we only consider relatively “red” blue GCs, say those 
with (g — i) o > 0.75, the resultant t; rms of the 22 GCs 
between 30' and 60' is 252+gg km s -1 , whereas for the 
other 32 blue GCs with (g — i) o < 0.75 the corresponding 
velocity dispersion is 524+34 km s _1 . Since lower lumi¬ 
nosity galaxies have on average bluer colors than higher 
luminosity galaxies for both their blue and red GC sys¬ 
tems (e.g. Peng et al. 2006), it is quite plausible that the 
high-dispersion “bluer” GCs are overwhelmingly contam¬ 
inated in projection by the intra-cluster populations that 
have been tidally stripped from dwarf galaxies, whereas 
the relatively “redder” GCs trace the underlying poten¬ 
tial of M87 more faithfully. This suggests that the stellar 
halo of M87 extends beyond, instead of being truncated 
at (Dopherty et al. 2009), R av ~ 150 kpc. 

7. VELOCITY DISTRIBUTION 

The line-of-sight velocity distribution is shaped by the 
global kinematics of a system. FigureHUlpresents the line- 
of-sight velocity histograms (with the systemic velocity 
of M87 being subtracted) of UCDs (left), blue GCs (mid¬ 


dle) and red GCs (right). We show the bright and faint 
subsamples in Figures [Tl] The bin size of the histograms 
is 50 km s -1 , which was chosen to be larger than the 
typical measurement uncertainties. The red curve over¬ 
plotted on each histogram is the adaptive kernel density 
estimation (KDE) of the velocity distribution. The KDE 
was constructed through a Gaussian-shaped kernel, with 
the Gaussian a being equal to the measurement uncer¬ 
tainty for each data point, and we further smoothed the 
resultant KDE with a Gaussian of a = 20 km s _1 (~ 
1/2.35x the bin size of the histograms). Also overplot¬ 
ted on each histogram as dashed blue curve is a scaled 
Gaussian distribution, with the Gaussian a being equal 
to the standard deviation of the observed distribution. 

To quantify the overall shape of the velocity distri¬ 
bution, we calculated the standard deviation a, skew¬ 
ness Gi, and the kurtosis G 2 , and they are listed on 
the top of each panel in Figures [101 and [CD The skew¬ 
ness was calculated as the ratio of the third cumulant 
and the 1.5th power of the second cumulant, and the 
kurtosis was calculated as the ratio of the fourth cumu¬ 
lant and the square of the second cumulant. For small 
samples, skewness and kurtosis defined by cumulants are 
relatively unbiased compared to the traditional defini¬ 
tions with moments. A Normal distribution has both the 
skewness and kurtosis equal to zero, and a distribution 
with sharper peak and (especially) heavier tails has more 
positive kurtosis. 

The estimation of kurtosis is sensitive to extreme out¬ 
liers. A meaningful estimation of kurtosis should reflect 
the overall shape of a distribution, instead of being driven 
by few extreme outliers. To obtain a robust estimation of 
kurtosis, we adopted a “cr-clipping” method in the kur¬ 
tosis space. Specifically, for a data set of N points, we 
resampled the original data set N times (without replace¬ 
ment), with one different data point being taken out each 
time, similar to the technique of jackknife resampling. If 
the resultant kurtosis after taking out a given data point 
is more than 5a away from the mean, that data point 
is regarded as an outlier. This “cr-clipping” process was 
iterated until no further outlier was found. Among the 
full samples within R < 30' (Section [TT]), 2 UCD (2%), 
13 blue GCs (2%) and 5 red GC (2%) were found to be 
outliers for kurtosis estimation. The clipped outliers are 
mostly the extreme velocities in our samples, and it is 
quite possible that most of these outliers belong to the 
intra-cluster population of the Virgo cluster. 

As a fourth-moment measurement, it is not surprising 
that the standard kurtosis largely reflects the tail be¬ 
havior. A complete description of the shape properties 
of a distribution should involve both the tailedness and 
peakedness. An outlier-insensitive, quantile-based alter¬ 
native for the standard kurtosis, i.e. the T parameter, 
was introduced by Moors (1998), and this alternative 
definition is expected to be more sensitive to the peaked¬ 
ness than Gi. A detailed introduction about T is given 
in the Appendix. The T parameter is defined such that 
a Normal distribution has a T equal to 0, and a positive 
T indicates heavier tails and (especially) a sharper peak 
than a Normal distribution. 

The a , G\ and G 2 reported below were calculated 
based on outlier-rejected samples. The 68% confidence 
intervals for all the above mentioned shape parameters, 
including the T parameter, were determined with by ran- 
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Fig. 7.— The line-of-sight velocities vs. the geometric average radius from M87. The UCDs, blue GCs, red GCs, and dE galaxies 
(non-nucleated and nucleated) are plotted separately as different symbols. The horizontal dotted line marks the systemic radial velocity of 
1307 km s -1 for M87 (Binggeli et al. 1993). 
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Fig. 8.— The UCD line-of-sight velocities vs. the geometric average radius from M87. Here different UCDs are color-coded with their io 
magnitude. The horizontal dotted line marks the systemic radial velocity of 1307 km s —1 for M87. 
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domly resampling the real data sets. The estimated G 2 
and T for some specified radial bins are listed in Tables 

[U ED and 01 


7.1. The Full Samples 

The UCDs and blue GCs have similarly higher dis¬ 
persion in the velocity distribution than the red GCs. 
The velocity distribution of UCDs is skewed toward the 
higher velocity tail, as quantified by a positive skewness, 
whereas the distribution of red GCs is skewed toward the 
lower velocity tail, as quantified by a negative skewness. 
The skewness difference between the three populations 
is partly reflected in their different systemic velocities 
T4 ys (Table [2]). The velocity distribution of UCDs is no¬ 
ticeably sharper than a Gaussian (blue dashed curves 
in Figure ITUll . whereas the velocity distributions of GCs 
are only marginally sharper than a Gaussian. The dif¬ 
ferent sharpness of the three distributions is well re¬ 
flected in their different T parameters. In addition, the 
UCDs and blue GCs have a similarly negative G 2 , sug¬ 
gesting slightly lighter tails than a Gaussian. As we will 
show later (Section [9]), the peaky velocity distribution of 
UCDs is consistent with a radially-biased velocity disper¬ 
sion tensor at large galactocentric distances, whereas the 
lighter tails are in line with a tangentially-biased velocity 
dispersion tensor at small distances. 

7.2. Bright and Faint Subsamples 

Since there are only 11 UCDs at io > 20.5 mag, the cal¬ 
culation of skewness and kurtosis for these faint UCDs 
is subject to large uncertainties and bias, and will not 
be discussed further. For the blue GCs, the bright and 
faint subsamples have similar velocity dispersion. Never¬ 
theless, the bright blue GCs have much more negative 
G 2 and marginally lower T than the faint ones. For the 
red GCs, the bright subsample has a significantly smaller 
(by ~ 70 km s -1 ) velocity dispersion and larger G 2 than 
the faint subsample. 

Most of the confirmed GCs with io > 21.5 mag were 
observed by Sll with the Low Resolution Spectrometer 
(LRIS) on Keck. The Keck/LRIS survey of Sll only cov¬ 
ered the central ~ 1.5 - 5.5' of M87, and as a result the 
confirmed faint GCs are primarily located in the central 
region (Figure [5j). To check if the velocity distribution 
differences between the bright and faint GCs are driven 
by the observational bias in spatial coverage for the sam¬ 
ple of faint GCs, we derived the shape parameters of the 
velocity distribution for the bright and faint GCs within 
the central 5'. It turns out similar differences between 
the bright and faint subsamples still exist for both the 
blue and red GCs. 

Furthermore, previous studies (Sll; Agnello et al. 
2014) found that the M87 GCs of different colors may ex¬ 
hibit different kinematical properties. So we also checked 
for any possible color bias for the bright and faint sub¬ 
samples. The median {g — i) 0 colors of the bright and 
faint blue GCs are 0.76 and 0.74 respectively, suggesting 
that there is no significant color bias for the bright and 
faint subsamples. For the red GCs, the median ( g — i) 0 
colors of the bright and faint subsamples are 0.95 and 
1.00 respectively. When we divide the red GCs into two 
(g — i)o color groups with a division color of 0.97, the 
above mentioned difference still exist for subsamples in 


each color group, although there is a systematic differ¬ 
ence between the “bluer” and “redder” groups, in the 
sense that the corresponding subsamples in the “bluer” 
group have about ~ 70 km s -1 lower velocity dispersion 
than those in the the “redder” group. 

7.3. Radial Trend of Kurtosis of the blue GCs 

The large sample size of the blue GCs allows us to 
explore the radial trend of the shape parameters for the 
velocity distribution. Figure [T2] shows the kurtosis for 
the bright and faint blue GCs in three fixed radial bins, 
namely 0'-4', 4'—12', and 12'-30'. It is significant that 
the bright subsamples have systematically lower kurtosis 
than the faint ones. 


7.4. Interpretation: from kurtosis to orbital anisotropies 

For a given number density profile, there is a one- 
to-one relationship between the kurtosis and the orbital 
anisotropy parameter /3 (= 1—crf/2of, where er* and oy 
are the tangential and radial components of the veloc¬ 
ity dispersion ellipsoid; positive (negative) /3 corresponds 
to radially (tangentially) anisotropic orbital structures), 
under the assumption that the velocity dispersion and 
/3 are constant with galactocentric radius (Napolitano et 
al. 2009). Napolitano et al. (2009) found that /3 can be 
directly estimated through deprojection of observables 
involving the kurtosis of V[ os and volume number density 
profiles of the tracer population (see the formulae B10, 
Bll and B12 in Napolitano et al. 2009). 

With the kurtosis and surface density profiles 0 (pa¬ 
rameterized as Sersic functions, Section |5~41) in hand, 
the average < /3 > for UCDs, blue GCs and red GCs 
were estimated (through numerical integration from one 
effective radius to infinity for each of the three popula¬ 
tions, also listed in Tables and mD to be —0. 3 + 0 ) 4 , 

—0.3to.2> an d respectively for the full samples. 

For the bright subsamples, the corresponding < ft > are 
—0. 3 ^ 0 ) 4 , — l.Otg^, and 0.21 q 2 for the UCDs, blue GCs 
and red GCs respectively; For the faint subsamples, the 
< /? > are estimated to be —O.IIq .2 an d 0.7^0.1 for the 
blue and red GCs respectively. We emphasize that the 
adopted approach to calculate /3 is only strictly applica¬ 
ble to a constant velocity dispersion profile which is not 
true for either of our three populations, so our estimation 
of /3 here should be regarded at most as a zeroth-order 
approximation of the average < j3 >. Later in this pa¬ 
per (Section [U]), we will solve the Jeans equations for the 
radial profiles of /3. 

Taken at face value, among the three populations, the 
red GCs have on average the highest radial anisotropy, 
whereas the blue GCs have the lowest radial anisotropy. 
In addition, the faint subsamples tend to have higher 
radial anisotropies than the bright subsamples. Given 
the good correspondence between kurtosis and < /? >, 
the radial kurtosis profiles shown in Figure [lU] indicate 


t To deproject the observed surface number density profiles N R 
to volume number density profiles n r , we assumed a spherically 
symmetric geometry, which gives 


n r oc 


/" 


dN R 


dR 


dR V7? 2 — r' 2 


(4) 
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Fig. 9. — Radial variations of the rms velocities of UCDs, blue GCs, red GCs and dE galaxies. The profiles were constructed by sliding 
bins of fixed radial widths (A R av = 4' for R aw < 8'; A R av = 8' for 8'< Raw < 20'; A R av = 12' for 20'< R aw < 40'; A R av = 16' for R aw > 
40 / ), with a step of 1' between adjacent sliding windows. The solid curves correspond to the calculated u rms , and the shaded regions mark 
the 68% confidence limits. To be used in constructing the profiles, at least 19 data points should be available for a given radius bin. 
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Fig. 10.— The line-of-sight velocity distribution of the UCDs (left), blue GCs (middle) and red GCs (right) in the inner 30' of M87. The 
histograms have a bin width of 50 km s —1 . The red solid curve overplotted on each histogram is the kernel density estimation of the observed 
distribution, and the blue dashed curve represents a Gaussian distribution, with the Gaussian <j being equal to the standard deviation of 
the velocity distribution. Note that the Gaussian distributions have been scaled by the area under the histograms. The standard deviation 
cr, skewness Gi, kurtosis G 2 and T parameter (see the Appendix for the definition) were also listed on the top of each panel. Note that, the 
UCDs have a sharper velocity distribution than the Gaussian, whereas the GCs (especially the blue GCs) have overall close-to-Gaussian 
distributions. 
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Fig. 11. — The line-of-sight velocity distribution of the subsamples of bright (lower panels) and faint (upper panels) UCDs (left), blue 
GCs (middle) and red GCs (right). The histograms have a bin width of 50 km s —1 . The red solid curve overplotted on each histogram is 
the kernel density estimation of the observed distribution, and the blue dashed curve represents a Gaussian distribution, with the Gaussian 
a being equal to the standard deviation of the velocity distribution. Note that the Gaussian distributions have been scaled by the area 
under the histograms. The standard deviation cr, skewness G i, kurtosis G 2 and T parameter were also listed on the top of each panel. 


i? av (kpc) 



Fig. 12.— Radial variation of velocity Kurtosis of the bright 
(large open circles) and faint (small filled circles) blue GCs. 


fit to given observations (e.g. Vi±Avi) is: 



1 (Vi l/mod) 
."2 + . ' 


(5) 


where a p is the intrinsic (projected) velocity dispersion, 
and 

'Cmod = Uys “l” ''-'VotSil 1 (9 $o) (6) 

with 0o being the position angle (PA, measured east of 
north) of the rotation axis. Maximization of Equation [4] 
is equivalent to minimizing the y 2 statistic (e.g. Bergond 
et al. 2006): 


X 


2 


\ ' f (Vi C mo d) 
i l CT p + ( Av i ) 2 


+ In kp + (Aci) 2 ] | 


(7) 


that the faint blue GCs are more radially-biased than the 
bright ones. 

8. KINEMATIC MODELING 
8.1. Methodology 

Under the assumption that the intrinsic angular ve¬ 
locity of a system, of either GCs or UCDs, is only a 
function of the galactocentric distance r, the (projected) 
average line-of-sight velocities vary sinusoidally with the 
projected azimuth 9 (e.g. Cote et al. 2001). As a conse¬ 
quence, it is customary to fit the line-of-sight velocities 
of GC and UCD systems with sine or cosine curves as a 
function of projected azimuth 9 , in order to determine 
the rotation amplitude and rotation axis (e.g. Cohen & 
Ryzhov 1997; Kissler-Patig & Gebhardt 1998; Cote et al. 
2001; Sll; Pota et al. 2013). Assuming Gaussian distri¬ 
butions for both the measurement uncertainties and the 
intrinsic velocity dispersions, the likelihood of a model 


By minimizing Equation [7] for a given dataset, we can 
determine the systemic velocity v sys , rotation amplitude 
v TO t , rotation axis 9q and the intrinsic velocity dispersion 
tip. Note that we have assumed the kinematic axis ratio q 
to be 1 when writing Equation[6]in order to be consistent 
with previous kinematical studies of M87 GCs (e.g. Cote 
et al. 2001; Sll). The results are essentially the same 
when fixing q to the photometric axis ratio of the diffuse 
stellar light. Considering the possible (unknown) incli¬ 
nation of the rotation axis with respect to the plane of 
sky, u rot determined here should be regarded as a lower 
limit. To estimate the uncertainties of the fitted param¬ 
eters, we randomly resampled the real dataset and re¬ 
peated the kinematic fitting to the resamples, and then 
obtained the 68% confidence intervals from the resultant 
parameter distribution. 

8.2. Bias and Significance 
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Kinematics fitting to discrete data points tends to over¬ 
estimate the intrinsic rotation, and the degree of overesti¬ 
mation or bias depends on the sample size, the azimuthal 
distribution of data points, and the importance of rota¬ 
tion as compared to the dispersion (e.g. Sharpies et al. 
1998; Romanowsky et al. 2009; Sll). Reliability of the 
fitted rotation amplitudes can be quantified in two mu¬ 
tually related ways, one is the most likely level of bias, 
and the other one is the significance (or confidence level, 
CL) of the fitted rotation. 

Bias. To check the level of bias for the best-fit rotation, 
we first constructed a series of kinematic models for each 
of the three populations, with the rotation amplitudes 
varying from 0 to 250 km s -1 and the other kinematic 
parameters being fixed at their best-fit values from the 
real data sets. Then, starting from each of these kine¬ 
matic models, velocities at each observed GC or UCD lo¬ 
cation were randomly drawn from a normal distribution, 
with the variance being equal to a quadrature combina¬ 
tion of the velocity dispersion and measurement uncer¬ 
tainties. In particular, at each input rotation amplitude, 
5000 Monte Carlo samples were generated for each pop¬ 
ulation, and the standard kinematic fitting (Section 18.11) 
was carried out for these mock samples. Lastly, for each 
input rotation amplitude, we determined the median of 
the corresponding 5000 best-fit output rotation ampli¬ 
tudes, and this defines a one-to-one relation between the 
input and the most likely output , which can be used to 
estimate the most probable bias in our best-fit rotation 
for the real data sets. 

Significance, we follow the procedure first introduced 
by Sharpies et al. (1998) to estimate the CL of our best-fit 
rotation amplitudes (see also Romanowsky et al. 2009). 
In particular, for each of the three populations, we ran¬ 
domly shuffled the position angles of the observed data 
points for 5000 times, and repeated the kinematics fitting 
to each realization. A random shuffling of the position an¬ 
gles can erase (if any) signatures of any possible rotation. 
Therefore, if the percentage of random realizations that 
lead to fitted rotation amplitudes greater than or equal 
to the best-fit value u ro t for the real data set is p, then 
the confidence level of u ro t can be estimated as 1 — p. In 
this paper, confidence levels > 90% are regarded to be 
significant. 

8.3. Global Kinematics 

8.3.1. The Full Samples 

Figure [Rj] shows the azimuthal variation of line-of-sight 
velocities for all M87 UCDs within the central 30', along 
with the blue GCs and red GCs for comparison. Overlaid 
on the data in each panel is the best-fit sine curves. The 
fitting results are summarized in Table [2] The UCDs and 
blue GCs have similar intrinsic velocity dispersion. The 
rotation amplitude of UCDs is more than 4 (2) times 
stronger than that of the blue (red) GCs. Additionally, 
the rotation axis of UCDs is roughly orthogonal to that of 
blue GCs. It is interesting that the red GCs, which have a 
smaller velocity dispersion than the UCDs and blue GCs, 
have a rotation axis that is more aligned with the UCDs 
rather than the blue GCs. Our best-fit parameters for 
the blue and red GCs are generally consistent with those 
determined by Sll (see their Table 14) within the mutual 
uncertainties. Note that the rotation angles reported by 


Sll are the direction of maximum rotation amplitude, 
which is 90° offset from the angular momentum vector. 
In addition, the kinematics parameters reported by Sll 
were already bias-corrected. 

Following the procedure described in Section 18.21 we 
estimated the most likely bias and the significance (or 
CL) of our best-fit rotation amplitudes, and the results 
are shown in Figure [14] and Table [2] According to bias 
test, the intrinsic rotation amplitude for UCDs is most 
probably overestimated by ~ 10 km s -1 , for blue GCs 
~ 6 km s -1 , and for red GCs ~ 1-2 km s -1 . As to the 
confidence levels, the probability that we found a rota¬ 
tion amplitude greater than or equal to the best-fit value 
for UCDs purely by chance is ~ 2% (CL = 98%), the 
probability is ~ 27% (CL = 73%) for blue GCs, and ~ 
12% (CL = 87%) for red GCs. 

Surface fitting to the spatial distribution (within the 
inner 30') of line-of-sight velocities is presented in Fig¬ 
ure as a color-coded background for each of the three 
populatons. The data points over plotted on the fitted 
surface are also color-coded according to their individual 
line-of-sight velocities. In order to bring out details of 
the fitted surface, individual data points with velocities 
>1550 km s -1 or <1050 km s -1 were not distinguished 
in colors from those with velocities = 1550 km s -1 or 
= 1050 km s -1 respectively. The surface fitting was car¬ 
ried out with the Kriging technique as implemented in 
R package fields. The Kriging technique has been re¬ 
cently used for exploring the velocity field of GC sys¬ 
tems (e.g. Foster et al. 2013) and the galactic spatial 
distribution of metallicities (e.g. Pastorello et al. 2014). 
Instead of simply applying an inverse distance weighting, 
Kriging takes into account both the spatial configuration 
and covariances of the dataset when assigning weights to 
neighboring data points for interpolation. In this work, 
the spatial covariance is assumed to be an exponential 
function of separation distance, and the smoothing pa¬ 
rameter A was fixed. The global rotation axes determined 
from our kinematics fitting were indicated as red arrows 
in Figure H5l and the arrow length is proportional to the 
best-fit rotation amplitude. 

Our global kinematics fitting is essentially driven by 
the clustering trend of V(os along the azimuthal direction, 
and primarily reflects the velocity field in the central re¬ 
gion where most data points are located. Accordingly, 
as shown in Figure 1151 there is generally a good match 
between the direction of arrows and the fitted Kriging 
velocity field in the central regions. The Kriging surface 
fitting is driven by the overall data configuration and spa¬ 
tial covariance, and in particular, by definition it is not 
influenced/biased by spatial clustering of data points. Of 
the three populations, the blue GCs have an overall ve¬ 
locity field that seems to be more closely aligned with the 
photometric major axis, while the UCDs have a velocity 
field in better agreement with that of the red rather than 
the blue GCs. 

The significantly stronger rotation of UCDs as com¬ 
pared to the GCs suggests that the UCDs are kinemat¬ 
ically distinct from the GCs. The smaller velocity dis¬ 
persion of red GCs is in agreement with the fact that 
red GCs are more centrally concentrated than the blue 
GCs and UCDs (Figure Similar to the blue GCs, a 
more or less minor axis rotation was also recently de¬ 
tected from integrated stellar-light spectra in the central 
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one effective radius of M87 (e.g. Arnold et al. 2013). In 
addition, a twist of the velocity field within the central 
half arc minute was recently reported by Emsellem, Kra- 
jnovic & Sarzi (2014) based on IFU spectra of the inte¬ 
grated stellar light. The misaligned velocity field across 
different radii and among different kinematic tracers (e.g. 
UCDs vs. GCs) all suggest that the halo of M87 is most 
probably triaxial, instead of being axisymmetric (e.g. 
Schwarzschild 1979; Statler 1991; Franx, Illingworth & 
de Zeeuw 1991; van den Bosch et al. 2008; Hoffman et 
al. 2009; Emsellem, Krajnovic & Sarzi 2014). 

8.3.2. The Bright and Faint Subsamples 

In this subsection, we explore the possible differences 
in global kinematics between the bright and faint sub¬ 
samples. The relevant kinematic fitting results are sum¬ 
marized in Tables [3] and 0] Figure [16] shows the azimuthal 
variation of line-of-sight velocities for the bright and faint 
subsamples of UCDs and GCs separately. Again, the 
best-fit sine curves are overlaid on the data points. For 
the blue GCs, no significant rotation was found for either 
the bright or faint ones. The bright red GCs exhibit a sig¬ 
nificantly lower (by ~ 50 km s -1 ) velocity dispersion and 
much less significant rotation than the faint subsample. 
The sample size of faint UCDs is too small to give mean¬ 
ingful kinematic parameters. The kinematic parameters 
of the bright UCDs resemble those of the full sample. 

8.4. Radial variation of Kinematics 

Radial variation of the kinematic parameters ($o, U-otj 
and er p ) is presented in FigurelT51 Fitting results for data 
points falling into some fixed radial range are listed in Ta¬ 
bles [3] and 01 The way that we constructed the profiles 
in Figure [18] is the same as in Section [6] which presents 
the w r ms profiles (Figure[9]). Basically, we performed kine¬ 
matic fitting to data points that fall into each individual 
sliding radial bin. Figure 07] shows the azimuthal varia¬ 
tions of line-of-sight velocities of UCDs (the left column), 
blue GCs (the middle column) and red GCs (the right 
column) in five different elliptical annuli. The five annuli 
were selected to be representative of the key features in 
the radial profiles of the UCD kinematics (Figure 05]). 
Note that the kinematic parameters plotted in Figure 05] 
are not corrected for possible bias. 

Based on Figure 05] and Table [2] the UCDs and blue 
GCs have similar er p profiles. er p of the red GCs is system¬ 
atically lower than the other two populations across the 
full radius range. The strongest rotation of UCDs is found 
around R av ~ 8'- 16', where the rotation axis is similar 
to the full sample of UCDs. The blue GCs exhibit in¬ 
significant rotation across the full radius range, whereas 
the red GCs show marginally significant rotation around 
~ 10'. We note that, in the Kriging maps shown in Figure 
01 a weak, but visible, gradient along the photometric 
major axis can be seen for blue GCs beyond the central 
~ 5', suggesting that the usual one-dimensional kinemat¬ 
ics fitting (Vj os vs. PA) and the two-dimensional surface 
fitting are complementary to each other. The seemingly 
strong rotation for blue GCs beyond the central 30' turns 
out to be not significant. 

The relevant kinematic parameters for the bright and 
faint subsamples of GCs in some specified radial bins are 
listed in Tables0]and0]separately. No significant rotation 


was found across the full radius range for both the bright 
and faint GC subsamples. 

Lastly, we checked the color-magnitude distribution 
of the UCDs in the radius range (10'-20') where the 
strongest rotation was found, and it turns out that these 
UCDs follow a color-magnitude relation indistinguish¬ 
able from the whole sample. This suggests that the strong 
rotation most probably reflects the global kinematics of 
the UCD system. 

8.5. Velocity Field of the Surrounding dEs 

The Virgo core region is still dynamically young (e.g. 
Binggeli et al. 1987). Deep wide-held optical imaging 
toward the core region revealed a complex network of 
extended tidal features surrounding M87 and other gi¬ 
ant ellipticals (e.g. Mihos et al. 2005; Janowiecki et al. 
2010), suggesting an ongoing hierarchical assembly of the 
Virgo core. It is thus natural to look for any connection 
between the velocity held of surrounding dwarf galaxies 
and stellar clusters in M87. Figure 09] presents the veloc¬ 
ity held for 69 non-nucleated (left panel), 59 nucleated 
(middle panel) and all (right panel) dE galaxies within 
2° of M87. As in Figure 05] Kriging maps of the mean 
velocity helds are color-coded in Figure 00] and the indi¬ 
vidual data points follow the same color scheme. There is 
no signihcant difference between the direction of velocity 
gradients of nucleated and non-nucleated dEs. The direc¬ 
tion of the velocity gradients of dEs more or less follows 
the photometric major axis of M87. Among the three ve¬ 
locity helds shown in Figure 05] the blue GCs seem to 
match the dEs best, in general agreement with Cote et 
al. (2001). 

We note that the remarkable velocity gradient on the 
Kriging maps of dEs is primarily driven by an excess 
of low-velocity dEs toward the north west direction. 
These excess low-velocity dEs are most probably associ¬ 
ated with a small subcluster of galaxies (e.g. Binggeli et 
al. 1993; Schindler, Binggeli, & Bohringer 1999; Jerjen, 
Binggeli, & Barazza 2004) centered on M86, which has a 
radial velocity ~ —244 km s _1 and is about 1 Mpc more 
distant than M87 (Mei et al. 2007). It is no doubt that 
the M87 subcluster and M86 subcluster are moving to¬ 
ward each other, and an imminent merging between the 
two of them has long been speculated (e.g. Bohringer et 
al. 1994; Binggeli et al. 1993). 

8.6. Interpretation: Ongoing Accretion of Dwarf 
Galaxies ? 

The observation that the surrounding dEs (especially 
the non-nucleated ones) follow a similar velocity field 
to the GCs (especially the blue ones) is consistent with 
the scenario that the GC systems, especially the metal- 
poor ones, may have been primarily assembled by ac¬ 
creting satellite dwarf galaxies along the photometric 
major axis of M87 (e.g. Cote, Marzke & West 1998). 
In line with this ongoing accretion or infalling picture, 
West & Blakeslee (2000) found that Virgo’s brightest 
ellipticals have a strong collinear arrangement in three 
dimensions. This so-called “principal axis”, which ap¬ 
pears to join a filamentary bridge of galaxies connecting 
the Virgo cluster to Abell 1367 and passes through the 
Virgo core, is also more or less aligned with the major 
axes of Virgo’s ellipticals (including M87). This “princi- 
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Fig. 13.— The line-of-sight velocity variations as a function of PA (East of North) for the full samples of UCDs (left), blue GCs (middle), 
and Red GCs (right). Overlaid on the data in each panel is the best-fit sine curve. The horizontal dotted lines mark the systemic radial 
velocity of 1307 km s —1 of M87, and the vertical dotted lines mark PA of the photometric minor axis of M87. 



Fig. 14. — Estimation of the most likely bias (left) and the confidence level (right) of the best-fit rotation amplitudes (Figures 1 13 1 and 
for UCDs, blue GCs, and red GCs at R :vv < 30'. In both the left and right panels, the black, blue and red colors correspond to results 
for the UCDs, blue GCs and red GCs respectively. In the left panel, curves of different colors correspond to the bias correction relation, as 
described in Section 18.21 for the three populations, and the open squares mark the best-fit rotation amplitudes for the real data sets. In 
the right panel, the open histograms represent the normalized distribution of best-fit rotation amplitudes for randomly shuffled samples, 
and the filled region of each histogram marks the probability p of obtaining a rotation amplitude greater than the best-fit value for the real 
data set just by chance, with the confidence level (CL) for the rotation being defined as 1 — p. 


pal axis” is thought to be the direction along which ma¬ 
terial flows into the cluster and forms galaxies, as seen 
in Cosmological IV-body simulations (e.g. van Haarlem, 
Frenk & White 1997; Hopkins, Bahcall & Bode 2005; 
Faltenbacher et al. 2005). 


9. ORBITAL ANISOTROPIES FROM JEANS ANALYSIS 
9.1. Method 


In this section we will infer the orbital anisotropies of 
UCDs and GCs based on the spherically symmetric Jeans 
equation 


GM(< r) 


d(n r a%) 

dr 



( 8 ) 


where M(<r) is the mass interior to the three- 
dimensional radius r. n r and ay are respectively the vol¬ 
ume number density and radial component of the veloc¬ 
ity dispersion at radius r for a given tracer population, 
and /3 r is the anisotropy parameter defined as 1—of/2of, 


with at being the tangential components (of = of + of) 
of the velocity dispersion (Binney & Tremaine 2008). 

In using Equation [SJ we have assumed that the net ro¬ 
tation can be either ignored or simply folded into the ve¬ 
locity dispersion term. This approximation is reasonable 
given the fact that the overall rotation of our samples is 
not dynamically important (t'rot/op < 0.4). In addition, 
as mentioned previously, the M87 system is probably tri- 
axial in shape, which would caution against a spherically 
symmetric Jeans analysis. However, as (at least) a first 
order approximation, it is definitely enlightening to do a 
comparative study of UCDs and GCs under the spheri¬ 
cally symmetric assumption. 

The anisotropy parameter /3 r can be constrained if 
M(< r), n r , and oy are known. For M(< r), we adopt 
the most recent determination by Zhu et al. (2014) based 
on made-to-measure modeling (Syer & Tremaine 1996; 
Long & Mao 2010) of over 900 M87 GCs, which extend 
out to a projected radius of ~ 180 kpc. The Zhu et al. 
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Fig. 15.— Spatial distribution of the UCDs {left), blue GCs ( middle ) and red GCs {right) are over plotted on their respective surface 
fitting {the color background) to the spatial distribution of line-of-sight velocities with the Kriging technique for the inner 30 / of M87. The 
black ellipses represent the stellar isophotes of M87 at 10R e , and the black solid line marks the photometric minor axis of M87 in each 
panel. The red arrows, with the length be ing proportional to the rotation amplitude, mark the direction of rotation axis from our global 
kinematics fitting to the inner 30' (Figure [13]) . The global kinematics fitting, which is primarily driven by the inner regions that contain 
most of the data points, matches the central velocity field from Kriging surface fitting. Among the three populations, the blue GCs seem 
to have an overall velocity field more aligned with the photometric major axis than the other two populations. See the text for details. 
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Fig. 16.— Radial velocity variations as a function of PA for the brighter (lower panels) and fainter (upper panels) subsamples of UCDs, 
blue GCs (middle), and Red GCs (right). Overlaid on the data in each panel is the best-fit sine curves. The horizontal dotted lines mark 
the systemic radial velocity of 1307 km s -1 of M87, and the vertical dotted lines mark PA of the photometric minor axis of M87. 


(2014) mass profile is a combination of a stellar compo¬ 
nent and a spherical logarithmic dark matter halo model. 
The n r profile is generally related to the surface number 
density profile Nr through the Abel integral equation. 
In particular, for deprojection of the Sersic surface pro¬ 
file that was used to characterize Nr of the UCDs and 
GCs in this work (Table [T]) , we adopt the analytical ap¬ 
proximation proposed by Prugniel & Simien (1997, Eq. 
B6; see also Mamon & Lokas 2005). In addition, treating 
the Jeans equation dHJ) as a first-order linear differential 
equation for n r a 2 , one finds (See also Cote et al. 2001; 
Mamon & Lokas 2005) 


n r e^ 2 ^ £*\l 


12 fS±dt] ni GM (< *) 


t 2 


dt (9) 


The line-of-sight velocity dispersion cr \ os ,r at projected 


radius R is determined (e.g. Binney & Mamon 1982) as 


c. 


o r°° d 2 

Jl _ ^ (- n 

los ,R 


= TT - fi r - T ) 


Nr 


r 2 y/r 2 - R 2 


dr (10) 


For the radial dependence of j3 r , we adopt the following 
function form (first proposed by Mamon, Biviano & Boue 
2013), 


fir fio T (fioo — fio)- 


rp 


( 11 ) 


where fio, fioo and rp are the three free parameters defin¬ 
ing the radial profile. Specifically, /? 0 and fioo are the 
anisotropies at r = 0 and oo respectively, and rp repre¬ 
sents the scale radius of fi r profile. This function form 
of fi r allows for either a radially increasing or decreasing 
profile. 

The GC (the blue plus red) anisotropy profile as de¬ 
termined by Zhu et al. (2014) shows a non-monotonic 
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Fig. 17.— Radial velocity variations as a function of PA in five elliptical annuli for the full samples of UCDs (left), blue GCs (middle), 
and red GCs (right). The dotted lines mark the systemic radial velocity of 1307 km s —1 for M87. The five elliptical annuli plotted here are 
selected to be representative of the primary features exhibited in the radial profiles of kinematics shown in Figure [18] 


behavior, in the sense that (3 r gradually increases to¬ 
ward the intermediate radii (~ 40 kpc) and then falls off 
in larger radii. Therefore, we also considered a two-part 
piecewise radial dependence of anisotropies by allowing 
the inner and outer radii to follow different profiles as 
defined by Equation II11 with the “transition” radius r tr 
being left as a free parameter. 

To constrain the radial anisotropy profiles for each of 
the different samples, we first created a library of model 
o"ios,R profiles for each of them by allowing the free pa¬ 
rameters that define the inner and outer anisotropy pro¬ 
files to uniformly (linearly for /3 q, /3oo, and rt r , logarith¬ 
mically for rp) vary. Then a maximum likelihood method 
was used to fit the models to the observed line-of-sight 
velocities as a function of projected radii. In particu¬ 
lar, by assuming that the observed line-of-sight veloci¬ 
ties Vi (±Ai>j) at a given projected radius R follow a 
Gaussian distribution with of = (An,) 2 + af os R and 
/r = we calculated a joint probability 

(similar to Eq.[5]in form) of each model profile for a given 
population. The most probable model profile is taken as 
the fiducial one, and the 68% confidence intervals are de¬ 
termined by randomly resampling the real data sets, with 
~ 10% of data points being left out for each resample. 


9.2. Results 

The derived anisotropy profiles for the full samples of 
UCDs and GCs within the central 35' of M87 are shown 
in Figure EHl Following Zhu et al. (2014), we reduced 
the weight of the data points (by increasing the uncer¬ 
tainties) that are located in the puzzlingly “hot” radius 
range from R av = 4' to 12' in the Jeans modeling. Since 
our Jeans analysis relies on the Zhu et al. mass profile, 
which was determined with made-to-measure modeling 
of nearly the same GC data set that is used in this work, 
we should obtain an anisotropy profile that is at least 
qualitatively consistent with Zhu et al. Comparing our 
Figure EU1 to the Figure 12 of Zhu et al. (2014), one can 
see that, although being based on different methods, our 
anisotropy profile for the blue GCs, which dominate the 
spectroscopic samples of M87 GCs, is in reasonably good 
agreement with Zhu et al. within the uncertainties. 

The UCD system has an anisotropy profile that be¬ 
comes more radial with radius, with f3 r being negative 
(tangentially-biased) within the inner ~ 20 - 40 kpc and 
being positive (radially-biased) beyond. We note that a 
radially-biased orbital structure for UCDs at larger radii 
is in line with a peaky velocity distribution shown in 
Section [3 The blue GC system has a radially increasing 
/ 3 r profile in the inner ~ 40 kpc but a radially decreas- 
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i? av (kpc) 



i? av (arcmin) 


Fig. 18.— Radial variation of the best-fit kinematic parameters. The profiles were constructed with sliding bins of fixed radial widths, 
as in Figure [9] The red (blue) dotted horizontal lines in the top panel mark the photometric semi-major (semi-minor) axes. To be used in 
constructing the profiles, at least 19 data points should be available for the kinematics fitting in a given radial bin. The solid curves are 
the best-fit kinematics parameters, and the shaded regions correspond to the 68% confidence limits. 


ing profile at larger radii. Among the three samples, the 
red GCs exhibit the largest radially-biased velocity dis¬ 
persion tensor across the explored radius range, which 
may be surprising but nevertheless in line with their rel¬ 
atively large velocity kurtosis (Figure [TOl Table [2]). We 
note that, although being based on different mass mod¬ 
els of M87, a highly radially-biased anisotropy (/? ~ 0.8 
at ~ 150 kpc) was also found in the outer stellar halo 
of M87 by Doherty et al. (2009) based on the integrated 
stellar absorption-line data at small radii and planetary 
nebulas (PNs; trace the stellar diffuse light, Coccato et 
al. 2009) kinematics at large radii (< 150 kpc), and this 
is in agreement with our finding for red GCs. @ 

§ Agnello et al. (2014) recently determined the anisotropies for 


354 GCs within ~ 100 kpc of M87, by dividing the GC system 
into three kinematically distinct subpopulations of different colors, 
i.e. blue, intermediate-color and red GCs, with intermediate GCs 
mostly being separated out from the classic red GCs. Agnello et 
al. found a mildly radially-biased anisotropy (~0.3) for their inter¬ 
mediate GCs and a slightly tangentially-biased anisotropy (~—0.2) 
for their red GCs. Nevertheless, a tripartition of M87 GCs may be 
still oversimplified. As was shown by Sll (c.f. their Figure 8), there 
is a complex color dependence of velocity dispersion of the classic 
red GCs, in the sense that the “intermediate” red GCs have a sig¬ 
nificantly higher dispersion than both the “bluer” and “redder” red 
GCs. So, there may be at least three kinematically distinct subpop¬ 
ulations for the classic red GCs alone, and the red GC system has 
not completely mixed dynamically. A complete understanding of 
the dynamics of the red GCs would have to wait for larger samples 
of radial velocities across the M87 system, in order to explore the 
full color dependence of their dynamics. 
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dE,non-Ns dE,Ns all dEs 



Aa(2000) (arcmin) Aa(2000) (arcmin) Aa(2000) (arcmin) 

Fig. 19.— Spatial distribution of the non-nucleated dEs (left), nucleated dEs (middle) and all dEs (right) are over plotted on their 
respective color-coded surface fitting to the line-of-sight velocities with the Kriging technique for the inner 120 / of M87. The black ellipses 
represent the stellar isophotes of M87 at 10 R e , and the black solid line marks the photometric minor axis of M87 in each panel. Note that 
the direction of velocity gradients for both the nucle ated and non-nucleated dEs roughly follows the photometric major axis of M87. When 
compared to the velocity fields shown in Figure [151 the dEs seem to be more aligned with the blue GCs than the UCDs and red GCs. 


TABLE 2 

Kinematics of the Full Samples 
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Note. — (1) The range of geometric galactocentric distances from the center of M87. (2) Number of data points used in the kinematics 
modeling. (3) The best-fit systematic velocity. (4) The best-fit azimuthal angle of the rotation axis, east of north. (5) The best-fit rotation 
amplitude. (6) The bias-corrected rotation amplitude. (7) The confidence level of the best-fit rotation. (8) The rotation-subtracted velocity 
dispersion. (9) The biweight root-mean-square velocity. (10) Standard kurtosis G 2 of the velocity distribution. (11) The T parameter 
proposed by Moors (1988). As a quantile-based alternative for the standard kurtosis, T is an increasing function of G 2 . See the Appendix 
for definition of T. 
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TABLE 3 

Kinematics of the Subsamples with io < 20.5 
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Note. — (1) The range of geometric galactocentric distances from the center of M87. (2) Number of data points used in the kinematics 
modeling. (3) The best-fit systematic velocity. (4) The best-fit azimuthal angle of the rotation axis, east of north. (5) The best-fit rotation 
amplitude. ( 6 ) The bias-corrected rotation amplitude. (7) The confidence level of the best-fit rotation. ( 8 ) The rotation-subtracted velocity 
dispersion. (9) The biweight root-mean-square velocity. (10) Standard kurtosis G 2 of the velocity distribution. (11) The T parameter 
proposed by Moors (1988). As a quantile-based alternative for the standard kurtosis, T is an increasing function of G 2 . See the Appendix 
for definition of T. 



r (kpc) 

Fig. 20. — Variation of the anisotropy parameters as a function of the 3D radius. The profiles for UCDs, blue GCs, and red GCs are 
represented as black, blue, and red solid curves respectively. Following the same color code, the hatched regions of different styles mark the 
68 % confidence intervals for blue GCs, and red GCs. The grey shaded region marks the 68 % confidence interval for the UCDs. The short 
dashed curves (black for UCDs, blue for blue GCs) represent the anisotropy profiles predicted by a universal relation between the number 
density slope and /3 for relic high-cr density peaks as found in cosmological simulations by Diemand et al. (2005). 
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TABLE 4 

Kinematics of the Subsamples with *o > 20.5 
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Note. — (1) The range of geometric galactocentric distances from the center of M87. (2) Number of data points used in the kinematics 
modeling. (3) The best-fit systematic velocity. (4) The best-fit azimuthal angle of the rotation axis, east of north. (5) The best-fit rotation 
amplitude. (6) The bias-corrected rotation amplitude. (7) The confidence level of the best-fit rotation. (8) The rotation-subtracted velocity 
dispersion. (9) The biweight root-mean-square velocity. (10) Standard kurtosis G 2 of the velocity distribution. (11) The T parameter 
proposed by Moors (1988). As a quantile-based alternative for the standard kurtosis, T is an increasing function of G 2 . See the Appendix 
for definition of T. 


9.3. Interpretation 

The different orbital anisotropies of the different pop¬ 
ulations may be attributed to either their different ori¬ 
gin or different orbital evolutionary histories. Our size- 
defined UCDs have at least an order of magnitude lower 
average density than GCs of similar luminosity, which 
means that UCDs are subject to a stronger tidal influ¬ 
ence when approaching small galactocentric distances. 
In addition, at a given average orbital radius, objects 
with more radially-biased orbits can plunge deeper into 
the central regions of their host, and thus are prone to 
stronger tidal disruption. Recent simulations (Pfeifer & 
Baumgardt 2013) suggest that UCD-sized clusters can 
be converted, via a continuous tidal stripping of their 
envelope, into GC-sized objects at small galactocentric 
distances. 

9.3.1. Depletion of Radial Orbits of UCDs at Small 
Galactocentric Distances 

The finding that UCDs in the central ~ 40 kpc are 
depleted of radial orbits may be partly attributed to 
a stronger tidal influence at smaller galactocentric dis¬ 
tances. To be quantitative, the Jacobi radius rj of a stel¬ 
lar cluster moving in the tidal field of its host galaxy is 
determined by its galactocentric distance r ga i, the galaxy 
mass Mg a i interior to r ga i and the cluster mass to. Specif¬ 
ically, rj = r ga i(m/2M gal ) 1 / 3 ( von Hoerner 1957; Inna- 
nen et al. 1983; Bertin & Varri 2008; Renaud et al. 2011; 
Webb et al. 2013). Based on this basic tidal theory and 
the Zhu et al. mass profile of M87, we estimated the ratio 



Fig. 21.— Radial variation of the ratio of limiting radius to 
Jacobi radius for M87 UCDs which have projected half-mass radii 
r h = 10 pc and M = 2x 1O 6 M0 (solid curves) or bxl0 G Mq (dashed 
curves). The thick and thin curves represent results respectively for 
r h/ r L — 0.07 and 0.14, which straddle the range of r^/r^ typical of 
the Galactic GCs based on King (1966) model fitting (McLaughlin 
Sz van der Marel 2005). The UCDs with \og(r^/rj) > 0 were tidally 
over-filling, and thus should have been tidally truncated below 10 
pc in rh- 

of limiting radius 7 *l to rj as a function of galactocen¬ 
tric distances for M87 UCDs (m = 2xl0 6 M 0 or 5xl0 6 
Mq) with projected half-mass radius r h = 10 pc (Figure 
ED- tl ideally represents the radius beyond which the 
stellar density of a cluster is zero, but in practice tt is 
an extrapolated radius from either King or other model 
profile fitting (e.g. Harris 1996). In plotting Figure I2T1 
we considered two different ratios of rh/rL, namely, 0.07 
and 0.14, which cover the typical range for the Galactic 
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GCs based on King (1966) model fitting (McLaughlin & 
van der Marel 2005). 

From Figure it can be seen that the limiting 
radii are comparable to or greater than the Jacobi radii 
(tidally over-filling) at galactocentric distances < 10 kpc 
for typical M87 UCDs if rh/n, < 0.07. Previous King 
profile fitting to some Virgo UCDs with the high-quality 
HST imaging data suggests that they have r^/r^ < 0.1 
(e.g. Evstigneeva et al. 2007). It is thus quite plausible 
that many UCDs with ?'h > 10 pc at small galactocen¬ 
tric distances have been tidally truncated into GC-sized 
objects with r h < 10 pc, and this tidal transformation 
should be especially efficient for UCDs on plunging orbits 
at a given average galactocentric distance. Accordingly, 
UCDs that survives at smaller galactocentric distances 
are expected to have less radially-biased orbital struc¬ 
ture, in line with our finding of an anisotropy profile that 
increases with radius for UCDs. 


9.3.2. Different Assemby Histories of UCD and blue GC 

Systems 

Although the tidal transformation effect discussed 
above should have been playing some role, it is not nec¬ 
essarily the sole or even dominant mechanism in shap¬ 
ing the present-day anisotropy profile of our size-defined 
UCDs. It is possible that the progenitors of many UCDs 
at smaller galactocentric distances have been primarily 
accreted from less radially-biased orbits. 

In a ACDM Universe, a more radially-biased orbital 
structure for the outer stellar halos of early type galax¬ 
ies seems to be a general outcome from a hierarchical 
structure formation (e.g. Dekel et al. 2005; Diemand 
et al. 2005; Hansen & Moore 2006; Abadi, Navarro & 
Steinmetz 2006). It has been conjectured that old halo 
populations, such as the metal-poor (blue) GCs, might 
form primarily in small dark matter halos that collapsed 
from high-cr (> 2.5cr; Diemand et al. 2005) peaks of the 
primordial density field at z > 10 (e.g. Moore et al. 
2006; Boley et al. 2009; Moran et al. 2014; see also a 
review by Brodie & Strader 2006). One interesting find¬ 
ing from the V-body cosmological simulation of Diemand 
et al. (2005) is the existence of a universal relation be¬ 
tween the present-day density profile slope dlnn/dhir 
and the anisotropy for the relic high-cr peaks, and this 
relation is not sensitive to the detailed assembly histo¬ 
ries. In particular, Diemand et al. (2005) found: /3 cs 
—0.23(1.2 + dlnn/dlnr). By deprojecting the surface 
number density profiles presented in Section 15.41 under 
the assumption of a spherically symmetric geometry, we 
obtained the /3 profiles for M87 UCDs and blue GCs as 
expected by the above relation. The profiles are overplot¬ 
ted in Figure [20] as dashed lines. 

Although being more negative in the innermost part, 
the fi r profile of UCDs determined from Jeans analy¬ 
sis is more or less consistent with the Diemand et al. 
predication within 2cr uncertainties. As discussed above, 
the significantly tangentially-biased orbital structure of 
UCDs at small radii can be partly attributed to a strong 
tidal transformation. The derived (3 r profile of the blue 
GCs exhibit large deviation from the Diemand et al. pre¬ 
diction. The finding that the blue GCs are tangentially- 
biased, rather than being radially-biased, at large radii 


may indicate^ that the blue GC system in the outer halo 
of M87 has not yet established an equilibrium state and 
is still in an early and active stage of assembly, pre¬ 
sumably through a continuous accretion of surrounding 
dwarf galaxies (e.g. Cote, Marzke & West 1998). Ag- 
nello et al. (2014) invoked the scenario of adiabatic con¬ 
traction of the dark matter halo to explain the enhance¬ 
ment of tangentially-biased orbits of blue GCs of M87 at 
large radii, and this mechanism works when accretion of 
satellite galaxies happens on sufficiently slow time scales 
(Goodman & Binney 1984). 

The agreement between the present-day radial 
anisotropy profile of UCDs and that predicted by cosmo¬ 
logical simulations does not tell us the detailed accretion 
process of the UCD progenitors. Our finding that UCDs 
have radially-biased orbital structure at large galacto¬ 
centric distances is in line with prediction of the “tidal 
threshing” scenario where UCDs were primarily tidally 
stripped dE nuclei (e.g. Bekki et al. 2003; Goerdt et al. 
2008; Pfeffer & Baumgardt 2013), although the primary 
progenitor galaxies of UCDs do not necessarily resem¬ 
ble the present-day surviving dEs. Previous simulations 
suggest that, to be tidally-threshed to a naked nucleus, 
a dE,N galaxy has to be on a highly radially-biased or¬ 
bit in order to plunge deep into the central potential. If 
tidal threshing is indeed the dominant channel for form¬ 
ing UCDs, our finding suggests that UCDs might be pref¬ 
erentially but not exclusively accreted from plunging or¬ 
bits. Indeed, previous simulations (e.g. Gill et al. 2004; 
Smith et al. 2013) suggest that, at small galactocentric 
distances, satellite galaxies on circular orbits may also be 
tidally disrupted. 

10. SUMMARY AND DISCUSSION 

We have compiled a sample of 97 spectroscopically con¬ 
firmed UCDs associated with the central cD galaxy of 
the Virgo cluster - M87 (NGC 4486). The UCDs are de¬ 
fined to have 10 < < 100 pc. 89% of our sample have 

io < 20.5 mag (M g < —10.6), corresponding to a stellar 
mass of ~ 2 x 1O 6 M 0 . In addition, 92% of the UCDs 
have colors as blue as the classic blue GCs. Throughout 
this paper, we compared the distribution and dynamics 
of UCDs to that of M87 GCs. The primary results of this 
paper are summarized as follows. 

1. The M87 UCD system has a surface number den¬ 
sity profile that is shallower than the blue GCs in 
the inner ~ 15' (~ 70 kpc) and as steep as the red 
GCs at larger radii. In addition, the UCDs follow 
a radial velocity dispersion profile more similar to 
that of the blue rather than the red GCs. 

2. Overall, the UCD system exhibits a significantly 
stronger rotation than the GC system. No sig¬ 
nificant rotation was found for either the bright 
(M g < —10.6) or faint (M g > —10.6) blue GCs. 
Although subject to relatively large uncertainties, 
the velocity field of dE galaxies surrounding M87 
seems to be more aligned with that of the blue GCs 
rather than the UCDs. 

^ Another possibility could be that the outer halo mass of M87 
may be somehow underestimated, which would artificially lead to 
a lower f3 due to the mass-anisotropy degeneracy. 
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3. The velocity distribution of UCDs is character¬ 
ized by a sharper peak and marginally lighter tails 
compared to a Gaussian distribution. This is in 
general agreement with results from Jeans analy¬ 
sis, namely, anisotropy of the velocity dispersion 
tensor of UCDs gradually increases from being 
tangentially-biased at inner radii to being radially- 
biased at larger radii. Overall, the GCs have veloc¬ 
ity distribution similar to a Gaussian, with the blue 
GCs being slightly platykurtic and the red GCs be¬ 
ing slightly leptokurtic. In addition, the bright blue 
GCs have velocity kurtosis systematically higher 
than the faint ones across the full range of galac- 
tocentric distances, indicating that the bright ones 
are more tangentially-biased than the faint ones. 

4. The M87 UCD system has an orbital anisotropy 
profile that gradually increases with galactocen- 
tric distances, with UCDs within ~ 40 kpc be¬ 
ing tangentially-biased while being radially-biased 
outward. In contrast to UCDs, the blue GCs be¬ 
yond ~ 40 kpc become gradually more tangentially- 
biased toward larger radii. The tangentially-biased 
anisotropy of UCDs in the inner region may be 
partly attributed to a continuous tidal transforma¬ 
tion of UCD-sized objects on plunging orbits to 
GC-sized clusters. 

Above all, our analysis suggests that the M87 UCDs 
are dynamically distinct from GCs. Other evidence 
against UCDs being the most luminous and extended 
tails of normal GCs include their different Fundamental- 
Plane relation (e.g. luminosity vs. internal velocity dis¬ 
persion, Cliilingarian et al. 2011). Our finding that UCDs 
have radially-biased orbital structure at large galacto- 
centric distances is in general agreement with the “tidal 
threshing” scenario that UCDs are primarily tidally 
stripped dE nuclei. In an accompanying paper, Liu et 
al. (2015) find that M87 UCDs at smaller galactocentric 
distances tend to have less prominent stellar envelopes 
than those lying at larger distances, providing direct evi¬ 
dence for tidal stripping. Moreover, previous studies (e.g. 
Paudel et al. 2010; Chilingarian et al. 2011; Francis et al. 
2012) suggests that Virgo UCDs have metallicities that 
are high for their luminosity according to the metallicity- 
luminosity relation defined by early-type galaxies, which 
is naturally expected for the “tidally threshed dwarf 
galaxy” scenario. 

In the context of ACDM hierarchical structure for¬ 
mation, as a more centrally confined population, UCDs 
might originate from rarer density peaks in the primor¬ 
dial density field than the more spatially distributed 
dwarf galaxies that have been presumably the main con¬ 
tributor of blue GCs to the outer halo of M87, and those 
rarer systems should collapse and fall into the central 
potential earlier. It is well known that the dE,N galax¬ 
ies are strongly centrally clustered in galaxy clusters as 
compared to the non-nucleated dEs (e.g. van den Bergh 


1986; Ferguson & Sandage 1989; Lisker et al. 2007). The 
distinct dynamical property of the UCD system might 
owe its origin to an earlier accretion of the progenitors 
of UCDs, and the present-day surviving dE,Ns do not 
necessarily resemble the primary UCD progenitors. In 
fact, most UCDs are found to be significantly older, more 
metal-poor and have super-solar alpha-element abun¬ 
dances compared to the majority of present-day dE nu¬ 
clei. Future spectroscopic stellar population analysis of 
large sample of M87 UCDs will be invaluable in further 
shedding light on the difference between the UCDs, blue 
GCs and dE nuclei. 
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APPENDIX 


A Quantile-based Alternative for Kurtosis 

As was detailed in the main text, velocity kurtosis is closely related to the orbital structure of a given tracer 
population. However, as a fourth moment measure, kurtosis in its standard form is more sensitive to heavy tails than 
to a sharper peak. A given distribution can be broadly divided into five parts, i.e. the peak, the two shoulders (e.g. /i 
± <r), and the two tails. The standard kurtosis measures the peakedness and tailedness of a distribution, which is more 
or less equivalent to measuring the dispersion (i.e. toward the peak and tails) around the shoulders, in the sense that 
a higher kurtosis indicates a larger dispersion around the two shoulders. Based on this interpretation, Moors (1988) 
proposed a robust quantile alternative to the standard kurtosis G 2 : 


T = 


(Et — E§) + (E 3 — Ei) 
Ee, — En 


- 1.23 


( 1 ) 


where £) is the i-tli octile. A normal distribution has T = 0. Moors (1988) showed that, although there is no simple 
relation between T and G 2 , T is an increasing function of G 2 . 

The T parameter, which is defined by quantiles, is more sensitive to peakedness than the standard kurtosis. The 
velocity distribution of our UCDs provides a good example to illustrate this point. UCDs have an obviously sharper 
peak than a Gaussian distribution, but this feature is only reflected in the high T value, not in the G 2 measurement. 


TABLE 5 

Ultra-compact Dwarfs 


ID 

R. A. (J2000) 

Decl. ( J2000) 

u los 

u* 

9 

r 

i 

z 

K s 

E(B- V) 

TijNGVS 

H^HST 


(deg) 

(deg) 

(km s- 1 ) 

(mag) 

(mag) 

(mag) 

(mag) 

(mag) 

(mag) 

(mag) 

(pc) 

(pc) 

(i) 

(2) 

(3) 

(4) 

(5) 

(6) 

(7) 

(8) 

(9) 

(10) 

(ii) 

(12) 

(13) 

Old Sample a 

H27916 

187.71521 

12.23610 

1299±10 

22.10:b0.01 

21.05±0.00 

20.52±0.00 

20.38±0.00 

20.22±0.01 

20.48±0.01 

0.024 

13.5±0.2 

13.7 

H30401 

187.82795 

12.26247 

1323±46 

22.54±0.01 

21.59±0.01 

21.06±0.00 

20.86±0.01 

20.74±0.01 

20.94±0.02 

0.022 

11.3±0.2 

10.7 

H30772 

187.74191 

12.26728 

1224±9 


20.75 


19.84 



0.023 


9.7 

H36612 

187.48603 

12.32538 

1601±3 

21.11±0.00 

19.99±0.00 

19.43±0.00 

19.16±0.00 

19.06±0.00 

19.16±0.00 

0.027 

17.5±0.3 

10.9 

H44905 

187.73785 

12.39440 

1563±18 

22.99±0.02 

21.92±0.01 

21.37±0.01 

21.14±0.01 

21.05±0.01 

21.13±0.02 

0.023 

22.3±1.1 

18.5 

H55930 

187.63929 

12.49845 

1297±4 

20.42±0.00 

19.24±0.00 

18.68±0.00 

18.45±0.00 

18.33±0.00 

18.41±0.00 

0.021 

32.9±0.5 

35.8 

S417 

187.75616 

12.32351 

1860±2 

21.00±0.00 

19.65±0.00 

19.03±0.00 

18.70±0.00 

18.49±0.00 

18.40±0.00 

0.023 

15.0±0.2 

14.7 

S477 

187.74961 

12.30030 

1651±62 

21.01±0.00 

20.06±0.00 

19.56±0.00 

19.34±0.00 

19.22±0.00 

19.52±0.01 

0.023 

23.8±0.5 

33.5 

S547 

187.73910 

12.42903 

714±2 

20.46±0.00 

18.87±0.00 

18.16±0.00 

17.76±0.00 

17.48±0.00 

17. lOibO.OO 

0.022 

20.3±2.5 

21.6 

S672 

187.72804 

12.36065 

735±106 

21.85±0.01 

20.83±0.00 

20.30±0.00 

20.08±0.00 

19.96±0.00 

20.18±0.01 

0.024 

19.3±3.0 

25.9 

S682 

187.72775 

12.33962 

1333±106 

22.21±0.01 

21.30±0.00 

20.81±0.00 

20.60±0.01 

20.50±0.01 

20.82±0.02 

0.024 

20.2±0.3 

23.7 

S686 

187.72421 

12.47187 

817±106 

21.53±0.00 

20.58±0.00 

20.05=b0.00 

19.83±0.00 

19.70±0.00 

19.97±0.01 

0.021 

16.3±0.6 

21.2 

S723 

187.72399 

12.33940 

1398±106 

22.68±0.01 

21.74±0.00 

21.21±0.00 

21.01±0.01 

20.85±0.01 

21.12±0.02 

0.024 


16.9 

S731 

187.72452 

12.28682 

1020±9 

22.24±0.01 

21. lOibO.OO 

20.56±0.00 

20.25ib0.00 

20.13±0.01 

20.19±0.01 

0.024 

20.7±0.4 

19.0 

S796 

187.71563 

12.34815 

1163±106 

21.84±0.01 

20.81±0.00 

20.28±0.00 

20.05±0.00 

19.92±0.00 

20.18±0.01 

0.024 

11.8±0.1 

15.3 

S825 

187.71263 

12.35542 

1142±106 

22.60±0.01 

21.63±0.00 

21.12±0.00 

20.91±0.01 

20.80±0.01 

21.17±0.02 

0.024 

12.8±1.0 

13.3 

S887 

187.70389 

12.36544 

1811±106 


21.19 


20.33 



0.024 


9.8 

S928 

187.69875 

12.40845 

1284±5 

20.81±0.00 

19.78±0.00 

19.26±0.00 

19.02±0.00 

18.90±0.00 

19.06±0.00 

0.023 

26.1±0.4 

23.0 

S999 

187.69130 

12.41709 

1467±5 

21.35±0.00 

20.30±0.00 

19.78±0.00 

19.52±0.00 

19.40±0.00 

19.58±0.00 

0.022 

20.6±0.6 

21.9 

S1201 

187.67423 

12.39478 

1211±106 

22.16±0.01 

21.18±0.00 

20.66±0.00 

20.42±0.00 

20.30±0.01 

20.60±0.01 

0.023 

14.5±0.4 

29.9 

S1508 

187.63087 

12.42356 

2419±140 

22.99±0.02 

22.01±0.01 

21.49±0.01 

21.31±0.01 

21.25±0.01 

21.53±0.03 

0.022 

22.2±0.4 

42.4 

S1629 

187.61066 

12.34572 

1129±7 

21.44±0.00 

20.38±0.00 

19.81±0.00 

19.63±0.00 

19.51±0.00 

19.69±0.01 

0.023 

18.0±0.3 

26.4 

S5065 

187.70854 

12.40248 

1578±3 

21.32±0.01 

20.21±0.00 

19.67±0.00 

19.40±0.00 

19.26±0.00 

19.33±0.00 

0.023 

12.4±0.3 

13.6 

S6004 

187.79259 

12.26697 

1818±77 

22.43±0.01 

21.32±0.00 

20.87±0.00 

20.56ib0.01 

20.43±0.01 

20.63±0.01 

0.022 

0.6±8.4 

40.3 

S8005 

187.69252 

12.40641 

1883±5 

21.56±0.01 

20.51±0.00 

19.97±0.00 

19.74±0.00 

19.61±0.00 

19.71±0.01 

0.022 

28.3±0.8 

25.9 

S8006 

187.69436 

12.40616 

1079±5 

21.61±0.01 

20.53±0.00 

19.99±0.00 

19.73±0.00 

19.62±0.00 

19.75±0.01 

0.023 

1.0±7.1 

21.2 

T15886 

188.15205 

12.34920 

1349±13 

23.96±0.03 

22.97±0.01 

22.55±0.01 

22.24±0.02 

22.22±0.03 

22.35±0.05 

0.028 

10.1±0.3 

11.0 

VUCD1 

187.53155 

12.60861 

1223±2 

20.28±0.00 

19.05±0.00 

18.51±0.00 

18.21±0.00 

18.04±0.00 

18.05±0.00 

0.022 

12.3=b0.1 

12.1 

VUCD2 

187.70085 

12.58636 

919±9 

20.29±0.00 

19.13±0.00 

18.57±0.00 

18.31±0.00 

18.17±0.00 

18.25ib0.00 

0.021 

11.1±0.1 

14.1 

VUCD4 

187.76865 

11.94347 

916±2 

20.30±0.00 

19.14±0.00 

18.64±0.00 

18.34±0.00 

18.20±0.00 

18.34±0.00 

0.028 

17.8±1.0 

25.1 

VUCD5 

187.79950 

12.68364 

1290±2 

20.44±0.00 

19.01±0.00 

18.38±0.00 

18.01±0.00 

17.82±0.00 

17.66±0.00 

0.025 

19.5±0.4 

19.2 

VUCD6 

187.86816 

12.41766 

2100±2 

20.47±0.00 

19.32±0.00 

18.76±0.00 

18.50ib0.00 

18.35±0.00 

18.45±0.00 

0.023 

13.1±0.3 

18.8 

VUCD7 

187.97040 

12.26641 

985±3 

19.76=b0.00 

18.49±0.00 

17.92±0.00 

17.58±0.00 

17.39±0.00 

17.38±0.00 

0.025 

19.6±0.4 

100.6 

VUCD9 

188.06074 

12.05149 

1323±12 

20.66±0.00 

19.45±0.00 

18.90±0.00 

18.57±0.00 

18.45=b0.00 

18.48±0.00 

0.029 

17.5±0.4 

25.4 


New Sample 


F6 

187.69749 

12.55047 

1341±5 

20.71±0.00 

19.46±0.00 

18.86±0.00 

18.56±0.00 

18.40±0.00 

18.34±0.00 

0.021 

17.4±0.2 

F12 

187.76079 

12.57058 

1190±14 

20.78±0.00 

19.79±0.00 

19.28±0.00 

19.05±0.00 

18.95±0.00 

19.18±0.01 

0.021 

18.2±0.2 

F16 

188.37279 

12.17151 

1230±18 

21.60±0.00 

20.57±0.00 

20.05±0.00 

19.81±0.00 

19.70±0.00 

19.91±0.01 

0.034 

12.2±0.3 
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TABLE 5 — Continued 


ID 

(1) 

R. A. ( J2000) 
(deg) 

(2) 

Decl. ( J2000) 
(deg) 

(3) 

”los 

(km s- 1 ) 

(4) 

(mag) 

(5) 

9 

(mag) 

(6) 

r 

(mag) 

(7) 

(mag) 

(8) 

(mag) 

(9) 

K s 

(mag) 

(10) 

E(B- V) 
(mag) 
(11) 

Hi,NGVS 
(pc) 
(12) 

Hi,HST 
(pc) 
(13) 

H18539 

187.51687 

12.12054 

1172±9 

20.84±0.00 

19.73±0.00 

19.20±0.00 

18.94±0.00 

18.85±0.00 

19.05±0.00 

0.028 

11. Iib0.2 


H20718 

187.58181 

12.15683 

861±9 

22.15±0.01 

21.08±0.00 

20.58±0.00 

20.38±0.00 

20.29±0.01 

20.58±0.01 

0.027 

11.0ib0.2 


H24581 

187.83334 

12.19947 

1283±9 

21.61±0.00 

20.48±0.00 

19.97±0.00 

19.71±0.00 

19.59±0.00 

19.74±0.00 

0.023 

11.4±0.2 


H51655 

187.94149 

12.45586 

1320±25 

22.30±0.01 

21.34±0.00 

20.86±0.00 

20.64±0.00 

20.52±0.01 

20.80±0.02 

0.025 

12.1±0.3 


H59533 

187.76536 

12.53695 

693±8 

21.04±0.00 

20.01±0.00 

19.51±0.00 

19.27±0.00 

19.16±0.00 

19.38±0.01 

0.020 

13.0±0.2 


H65115 

187.48801 

12.60081 

1491±10 

21.27±0.00 

20.27±0.00 

19.79±0.00 

19.54±0.00 

19.40±0.00 

19.69±0.01 

0.022 

12.8±0.2 


S41 

187.81801 

12.31261 

1790±31 

22.03±0.01 

20.94±0.00 

20.41±0.00 

20.14±0.00 

20.02±0.00 

20.14±0.01 

0.022 

26.2±0.4 


S323 

187.76845 

12.38923 

1157±106 

21.65±0.01 

20.67±0.00 

20.19±0.00 

19.95ib0.00 

19.83±0.00 

20.09±0.01 

0.023 

33.7±0.6 


S376 

187.75891 

12.46381 

1215±106 

21.59±0.01 

20.65ib0.00 

20.15±0.00 

19.92±0.00 

19.82±0.00 

20.09±0.01 

0.021 

23.4±0.3 


S804 

187.71277 

12.43663 

1137±7 

20.95±0.00 

19.67±0.00 

19.07±0.00 

18.75±0.00 

18.56±0.00 

18.50±0.00 

0.022 

13.5±0.3 


S991 

187.69376 

12.33826 

1004±75 

21.82±0.01 

20.79±0.00 

20.25±0.00 

20.00±0.00 

19.86±0.00 

19.98±0.01 

0.024 

15.5±0.3 


S1044 

187.68891 

12.34263 

2023±75 

21.70±0.01 

20.68±0.00 

20.14±0.00 

19.91±0.00 

19.78±0.00 

19.99±0.01 

0.024 

24.1±0.4 


S1301 

187.66316 

12.35901 

1086±106 

21.75±0.01 

20.72±0.00 

20.11±0.00 

19.91±0.00 

19.74±0.00 

19.96±0.01 

0.023 

13.6±0.4 


S1449 

187.64218 

12.37956 

1100±106 

22.37±0.01 

21.33±0.00 

20.80±0.00 

20.59±0.00 

20.46±0.01 

20.75±0.01 

0.022 

17.2±0.3 


S1504 

187.63137 

12.43405 

858±33 

21.26±0.00 

20.30±0.00 

19.75±0.00 

19.57±0.00 

19.45±0.00 

19.75±0.01 

0.022 

14.7±0.3 


S1617 

187.61207 

12.39196 

1407±28 

21.06±0.00 

19.93±0.00 

19.35±0.00 

19.13±0.00 

18.98±0.00 

19.08±0.00 

0.022 

16.7±0.2 


S1631 

187.60730 

12.43919 

1368±75 

21.52±0.00 

20.46±0.00 

19.93±0.00 

19.71±0.00 

19.59±0.00 

19.72±0.01 

0.022 

13.5±0.2 


S6003 

187.79226 

12.27445 

1818±77 

22.30±0.01 

21.21±0.00 

20.68±0.00 

20.40±0.00 

20.24±0.01 

20.42±0.01 

0.022 

20.3±0.4 


S9053 

187.70126 

12.49469 

829± 106 

22.53±0.01 

21.39±0.00 

20.89=b0.00 

20.60±0.00 

20.56=b0.01 

20.92±0.04 

0.021 

30.2±0.6 


VUCD8 

188.01813 

12.34176 

1647±3 

20.76±0.00 

19.64±0.00 

19.09±0.00 

18.82±0.00 

18.71±0.00 

18.82±0.00 

0.026 

12.4±0.3 


VUCD10 

187.62858 

12.31157 

2305±23 

20.86±0.00 

19.73±0.00 

19.18±0.00 

18.94±0.00 

18.79±0.00 

18.90±0.00 

0.023 

15.4±0.4 


M87UCD-1 

187.83029 

12.37554 

1136±21 

21.48±0.00 

20.34±0.00 

19.79±0.00 

19.53±0.00 

19.41±0.00 

19.51±0.00 

0.023 

28.4±0.3 


M87UCD-2 

187.69858 

12.14034 

1288±14 

20.77±0.00 

19.73±0.00 

19.25±0.00 

19.06±0.00 

18.90±0.00 

19.18±0.00 

0.026 

23.6±0.4 


M87UCD-3 

187.58354 

11.92186 

1404±13 

20.44±0.00 

19.34±0.00 

18.77ib0.00 

18.54±0.00 

18.37±0.00 

18.53±0.00 

0.030 

22.9±0.1 


M87UCD-4 

187.59096 

12.40067 

1279±11 

21.33±0.00 

19.93ib0.00 

19.29±0.00 

18.96±0.00 

18.75±0.00 

18.62±0.00 

0.023 

11.7±0.4 


M87UCD-5 

187.41829 

12.45780 

1400±27 

21.08±0.00 

19.96±0.00 

19.42:b0.00 

19.18±0.00 

19.08±0.00 

19.30±0.00 

0.023 

21.0±0.2 


M87UCD-6 

187.51025 

12.60938 

1514±27 

21.34±0.00 

20.21±0.00 

19.69±0.00 

19.42±0.00 

19.27±0.00 

19.37±0.01 

0.022 

15.0±0.1 


M87UCD-7 

187.52525 

12.66499 

1160±12 

20.65±0.00 

19.48±0.00 

18.95ib0.00 

18.67±0.00 

18.51±0.00 

18.58±0.00 

0.021 

17.4±0.1 


M87UCD-8 

187.43858 

12.85025 

1154±18 

21.60±0.00 

20.51±0.00 

19.99±0.00 

19.75±0.00 

19.63±0.00 

19.70±0.01 

0.020 

12.1±0.1 


M87UCD-9 

187.59571 

12.53021 

556±12 

20.87±0.00 

19.75±0.00 

19.22±0.00 

18.99±0.00 

18.89±0.00 

19.06±0.01 

0.022 

17.9±0.1 


M87UCD-10 

187.50812 

12.70747 

1178±30 

20.82±0.00 

19.86±0.00 

19.41±0.00 

19.17±0.00 

19.12±0.00 

19.35±0.00 

0.022 

21.9±0.3 


M87UCD-11 

187.62463 

12.63949 

845±30 

21.42±0.00 

20.47±0.00 

20.00±0.00 

19.79±0.00 

19.72±0.00 

20.02±0.02 

0.021 

19.0±0.6 


M87UCD-12 

187.63100 

12.86572 

1497±10 

21.51±0.00 

20.13±0.00 

19.48±0.00 

19.17±0.00 

19.01±0.00 

18.85±0.00 

0.021 

12.4±0.5 


M87UCD-13 

187.70533 

12.62171 

1605±7 

21.07±0.00 

19.95±0.00 

19.38±0.00 

19.12±0.00 

18.98±0.00 

19.01±0.01 

0.021 

12.8±0.2 


M87UCD-14 

187.76812 

13.17849 

1347±7 

21.64±0.01 

20.26±0.00 

19.67±0.00 

19.31±0.00 

19.13±0.00 

19.03±0.01 

0.025 

11.2±0.3 


M87UCD-15 

187.99162 

13.25948 

1261±12 

20.71±0.00 

19.66±0.00 

19.15±0.00 

18.90±0.00 

18.83±0.00 

19.08±0.01 

0.027 

40.2±0.4 


M87UCD-16 

187.88208 

12.69177 

1119±12 

20.65±0.00 

19.49±0.00 

18.89±0.00 

18.62±0.00 

18.46±0.00 

18.62±0.00 

0.025 

38.1±0.2 


M87UCD-17 

187.89638 

12.58316 

950±7 

21.40±0.00 

20. lOibO.OO 

19.47±0.00 

19.19±0.00 

19.00±0.00 

19.04±0.00 

0.024 

11.9±0.3 


M87UCD-18 

187.85313 

12.42384 

1780±8 

20.85±0.00 

19.80±0.00 

19.27±0.00 

19.04±0.00 

18.90±0.00 

19.09±0.00 

0.022 

13.2±0.1 


M87UCD-19 

187.99000 

12.82488 

1086±15 

21.38±0.00 

20.31±0.00 

19.81±0.00 

19.55±0.00 

19.48±0.00 

19.71±0.01 

0.023 

18.9±0.2 


M87UCD-20 

187.42217 

12.66457 

1754±105 

21.43±0.00 

20.25±0.00 

19.73±0.00 

19.45±0.00 

19.30±0.00 

19.42±0.01 

0.021 

16.7±0.3 


M87UCD-21 

187.90599 

12.29465 

1484±28 

21.59±0.00 

20.60±0.00 

20.16±0.00 

19.86±0.00 

19.77±0.00 

20.02±0.01 

0.023 

11.4±0.3 


M87UCD-22 

187.46783 

12.62716 

905±20 

22.09±0.01 

20.80±0.00 

20.21±0.00 

19.92±0.00 

19.72±0.00 

19.72±0.01 

0.021 

11.7±0.2 


M87UCD-23 

187.42317 

12.74100 

1314±23 

22.13±0.01 

20.95±0.00 

20.38±0.00 

20. lOibO.OO 

19.95±0.00 

20.04±0.01 

0.022 

11.1±0.2 


M87UCD-24 

187.29442 

12.71781 

1224±34 

22.01±0.01 

20.88±0.00 

20.32±0.00 

20.07±0.00 

19.97±0.01 

20.01±0.01 

0.022 

12.3±0.2 


M87UCD-25 

187.91054 

11.98996 

1264±26 

21.88±0.01 

20.78±0.00 

20.26±0.00 

20.01±0.00 

19.89±0.00 

20.13±0.01 

0.029 

19.4±0.4 


M87UCD-26 

187.56892 

12.26579 

2030±18 

21.05±0.00 

20.04±0.00 

19.55±0.00 

19.30±0.00 

19.17±0.00 

19.39±0.00 

0.023 

17.3=b0.3 


M87UCD-27 

188.10566 

12.33246 

1272±23 

21.56±0.00 

20.52±0.00 

20.01±0.00 

19.76±0.00 

19.69±0.00 

19.89ib0.01 

0.027 

11.7±0.3 


M87UCD-28 

187.52553 

12.40982 

1870±29 

21.78±0.00 

20.59ib0.00 

20.03±0.00 

19.74±0.00 

19.62±0.00 

19.68±0.01 

0.026 

12.9±0.2 


M87UCD-29 

187.02759 

12.41012 

599±33 

21.97±0.01 

20.94±0.00 

20.42±0.00 

20.21±0.00 

20. llibO.Ol 

20.32±0.01 

0.028 

12.6±0.3 


M87UCD-30 

187.15505 

12.47934 

1534±28 

22.06±0.01 

21.01±0.00 

20.48±0.00 

20.25±0.00 

20.14±0.01 

20.31±0.02 

0.024 

12.6±0.3 


M87UCD-31 

187.73054 

12.41109 

1301±30 

22.03±0.01 

20.95±0.00 

20.40±0.00 

20.16±0.00 

20.01±0.01 

20.08±0.01 

0.023 

9.7±0.1 

10.4 

M87UCD-32 

187.85521 

12.32549 

1632±34 

21.21±0.00 

20.22±0.00 

19.73±0.00 

19.48±0.00 

19.39±0.00 

19.65±0.00 

0.023 

33.2±0.5 


M87UCD-33 

188.06020 

12.03307 

1833±22 

21.81±0.01 

20.65±0.00 

20.09±0.00 

19.79±0.00 

19.66±0.00 

19.70±0.01 

0.029 

16.3±0.4 


M87UCD-34 

187.31196 

11.89551 

905±19 

21.25±0.00 

20.19±0.00 

19.65±0.00 

19.40±0.00 

19.27±0.00 

19.45±0.00 

0.027 

11.7±0.3 


M87UCD-35 

188.09787 

11.96527 

1007±16 

21.36±0.00 

20.19±0.00 

19.57±0.00 

19.34±0.00 

19.19±0.00 

19.23±0.00 

0.032 

11.5±0.3 


M87UCD-36 

187.79709 

12.50030 

1207±21 

22.03±0.01 

20.76ib0.00 

20.18±0.00 

19.88±0.00 

19.74±0.00 

19.79ib0.01 

0.020 

12.1±0.4 


M87UCD-37 

187.54091 

12.62679 

1324±38 

22.15±0.01 

21.04±0.00 

20.52±0.00 

20.26±0.00 

20.12±0.01 

20.20±0.01 

0.022 

11.3±0.1 


M87UCD-38 

187.62717 

12.67106 

1154±18 

21.31±0.00 

20.07±0.00 

19.49±0.00 

19.20±0.00 

19.09±0.00 

19.12±0.01 

0.022 

22.1±0.3 


M87UCD-39 

188.44296 

11.95064 

1351±23 

22.32±0.01 

21.00±0.00 

20.38±0.00 

20.03±0.00 

19.85±0.01 

19.77±0.01 

0.036 

11.2±0.2 



Note. — Col.(l): Object ID; Col.(2): Right ascension in decimal degrees (J2000); Col.(3): Declination in decimal degrees (J2000); Col.(4): Heliocentric radial velocity; 
Cols.(5—9): MegaCam u* griz five-band 3 // -aperture (in diameter) AB magnitudes (not corrected for Galactic extinction); Col.(10): WIRCam K s -band 3 // -aperture (in 
diameter) AB magnitude (not corrected for Galactic extinction); Col.(11): The Galactic reddening determined by Schlegel et al. (1998); Col.(12): Half-light radius (in 
units of pc) measured on NGVS images; Col.(13): Half-light radius (in units of pc) measured on HST images. 

The UCDs have been compiled by Brodie et al. (2011). In addition, g and i band photometry for H30772 and S887 is from Brodie et al. (2011) 






